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Poland 

We address the question of description of qubit system in a formalism based on the nilpo- 
tent commuting rj variables. In this formalism qubits exhibit properties of composite ob- 
jects being subject of the Pauli exclusion principle, but otherwise behaving boson-like. 
They are not fundamental particles. In such an approach the classical limit yields the nilpo- 
tent mechanics. 

Using the space of 77-wavefunctions, generalized Schrodinger equation etc. we study 
properties of pure qubit systems and also properties of some composed, hybrid models: 
fermion-qubit, boson-qubit. The fermion-qubit system can be truly supersymmetric, with 
both SUSY partners having identical spectra. It is new and very interesting that SUSY 
transformations relate here only nilpotent object. The 77-eigenfunctions for the qubit-qubit 
system give the set of Bloch vectors as a natural basis. 

Then the ?y-formalism is applied to the description of the pure state entanglement. 
Nilpotent commuting variables were firstly used in this context in [A. Mandilara, et. al., 
Phys. Rev. A 74, 022331 (2006)], we generalize and extend approach presented there. 
Our main tool for study the entanglement or separability of states are Wronskians of rj- 
functions. The known invariants and entanglement monotones for systems of n = 2, 3, 4 
qubits are expressed in terms of the Wronskians. This approach gives criteria for sepa- 
rability of states and insight into the flavor of entanglement of the system and simplifies 
description. 

Keywords: nilpotent commuting variables, nilpotent mechanics, entanglement, entanglement 
monotones, Stirling numbers, Bell numbers, symmetric polynomials, supersymmetic quantum 
mechanics, qubit-fermion systems 
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1. Introduction 

Qubit systems [1] are conventionally described in symmetric tensor products of the two 
dimensional complex Hilbert spaces. In such an approach there is no classical limit for 
the qubit understood as a particle. On the other hand, from the supersymetry theory we 
know that one can consider classical limit for fermions, at cost of introducing the new an- 
ticommuting variables. It turns out that similar goal can be achieved for qubits, but this 
time we have to introduce the nilpotent commuting variables. Such a new formalism is 
complementary to the conventional one, and gives natural setting to answer the entangle- 
ment questions. The aim of this work is to develop relevant formalism and apply it to the 
description of entanglement. 

As we already mentioned, one-particle space of the qubit is a two dimensional Hilbert 
space and its collective behavior is a boson-like. In terms of the Fock space it means that 
tensor product of qubit states is symmetric. Qubit from this point of view exhibits mixture 



Nilpotent quantum mechanics, qubits, and flavors of entanglement 3 



Table 1. Generalized commutation relations. 



Type 


Single site rel. 


Different site rel. 


Symmetry 


Bosons 
(CCR algebra) 


[b,b+]- = 1 
[b,b]- = 
[b+,b+]- =0 


[bi,b+]-=0 
[bi,bj]- =0 
=0 


+ 


Fermions 
(CAR algebra) 


[/,/+]+ = ! 
[/,/]+ =o 
[/+,/+]+ =o 


[/;,/+]+ =0 

[/*,/#]+ = o 

[/+,/+]+ =0 




Parabosons 


[e,e+]_ = 1 - 2\K 
[e,e]_ =0 
[e+,e+]_ =0 


[ei,et]_ =0 
[e«,ejj- = 
[e+,e+]_=0 


+ 


Parafermions 


[o,o + ]+ = 1 


[oi,ot]_ =0 


+ 


(Spin algebra) 


[o,o]+ = 
[a+,a+]+ =0 


[oi,Oj]_ = 
[a+ a+]_=0 




Qubits 

(Qubit algebra) 


= 1 - 2N 
d 2 = 
(d+) 2 = 


=0 

[di,dj]- = 
[d+ dt]-=0 


+ 



of the boson and fermion properties. Only separated qubit is like a fermion. To adequately 
describe qubits we have to play with the commutation or anticommutation relations on 
the one side, and symmetry properties of the tensor product on the other side. Bosons and 
fermions can be organized in unique graded structure where commutators, anticommuta- 
tors, parity of elements of graded algebra, and symmetry of tensor product is consistent. 
They play distinguished role, because they describe fundamental particles. However, there 
are other useful objects related to parastatistics. Parafermions and parabosons were defined 
by tri-linear relations in ref. [2] over fifty years ago and have their place in quantum field 
theory [3]. Here in our introductory discussion we use the reduced bi-linear form of these 
relations, that is satisfactory for our purposes, but less general then tree-linear relations. 
Such reduced form is used in the entries of the Table 1 concerning parafermions and para- 
bosons. A nice discussion of the definition of the qubit, stressing that qubits are neither 
bosons nor fermions, is given by Wu and Lidar [4]. We departure from their definition of 
qubit as a parafermion, but find it more useful to reformulate defining conditions using the 
commutator and nilpotency conditions. We collect defining properties of relevant objects in 
the Table 1 to compare various approaches with canonical commutation relations (CCR), 
canonical anticommutation relations (CAR), and parabosonic or parafermionic relations. 
In the last column there are given eigenvalues ±1 of the transposition (flip) map r, where 
T ii+ i(... ® ipi ® ip i+1 ...) = ... ® tpi+i <8> il>%— ■ The two dimensional state space of a sin- 
gle qubit system is naturally obtained from nilpotent creation/anihilation operators. Taking 
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commutators (symmetric tensor products) for many-site system we obtain parafermions. 
This approach is adopted most frequently cf. Ref.[4]. To find another description of qubits 
we rewrite the set of (anti)commutation relations for parafermions in the following way 

[d,d+]_ = l-2N, [di,dt]_=0, [di,dj]-=0, [dt,d+]_=ft.i) 

d 2 = 0, (d+) 2 = 0, 

where N can be seen as particle number operator. As it is known there is no classical limit, 
in the usual sense, for fermions. Only introduction of the Grassmannian variables allows 
to define such a limit. Analogously for qubits, to get nontrivial classical limit we have to 
introduce nilpotent, but this time - commuting variables. Namely, 

rp/ = r/v, v 2 = V /2 = (1.2) 

This provides that single qubit is two-level one, in the same time it is boson-like when 
considered in multi-qubit system. In the approach, where the qubits are realized within the 
conventionally defined spin algebra (parafermions), natural classical limit of such system 
yields the anticommuting co-ordinates [4], like for fermions. Such anticommuting coor- 
dinates imply naturally the antisymmetry of the tensor product, and nilpotent "bosonic" 
objects again have to be represented by even product of anticommuting variables. That is 
why parafermions in such a setting do not suit our demands. To describe qubit without 
referring to its composite character one has to use nilpotent commuting variables. But we 
have to pay for using the nilpotent ^-variables; the ^-"derivative" does not fulfil the Leib- 
niz rule [5, 6]. On the other hand, the system of many qubits in the //-formalism has the 
same property as observed in [4] for parafermions, that it behaves for large n like boson, 
in paricular becames not nilpotent. Already in the 1994 Palumbo [7] noted, in the context 
of high energy physics, that nilpotent commuting variables: 

[nilpotent commuting variables] are not just abstract mathematical entities devoid 
of any physical interpretation, but describe composites of fermions, so that models 
can be of phenomenological importance. 

The formalism of //-functions which later will be used to describe the nilpotent quan- 
tum mechanics was discussed in more detail in [5, 8, 9]. Properties of such functions allow 
to answer many questions related to the entanglement. Criteria of factorization of states 
coming from this formalism are natural, are derived by simple argument and are equiva- 
lent to the ones known from the invariants theory approach. In the work of Wu and Li- 
dar there is analysis of the notion of qubit from the particle point of view and is given 
the parafermionic description of this object. In our approach we rewrite the parafermionic 
character of qubit in such a way that the commutators and nilpotency conditions are used in 
defining relations, but resulting properties are analogous as in [4]. Despite presented above 
fundamental arguments for representing qubits with use of the ^-variables, arguments re- 
lated to symmetry of qubit state product and nilpotency; it was found by Mandilara et al. 
in [10, 11] that nilpotent commuting variables are very useful tool in the description of 
entanglement. One may think that such formalism has deeper roots, and like for fermions 
one can consider (pseudo)classical systems of superparticles, then quantize them. Here we 
should have the formalism of nilpotent mechanics for systems that after quantization will 
give the (quantum) qubit systems. 

Nilpotent commuting variables, besides the quantum mechanics [10], are present in 
the theoretical physics in several contexts: nuclear and high energy physics, string theory 
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[12, 13, 14, 7, 15], as well as in metioned before classical nilpotent mechanics. For review 
of some applications cf. Ref. [9]. 

The paper is organized as follows. In the next section we describe the qubit system in 
such a way that it has the proposed above classical limit realized by the commuting nilpo- 
tent coordinates rj. Then necessary formalism based on 77-coordinates is presented. We re- 
call definition of the 77-numbers, 77-functions [5] and then develop it by studying properties 
of the elementary ^-functions, in particular: exponent, logarithm, and trigonometric func- 
tions. It turns out that the Stirling numbers, the ordered Bell numbers, and the second kind 
Bell numbers are related to the expansions of specific ry-functions. The following section 
we devote to the introduction of the symmetric ^-polynomials and Hermite ^-polynomials. 
Then the rudiments of the f7-calculus are given with the ^-integral, the ^-Fourier trans- 
form as well as the ^-integral form of the Stirling and the normal Bell numbers. Having 
above tools we introduce the notion of the A^-Hilbert space which will allow to describe 
qubit systems. Using the //-integral we define the A/"-scalar product in the space of the 
ry-functions (wave 77-functions) and the structure of the 7V-Hilbert space. It turns out that 
natural bases in such spaces are closely related to the ones used in the description of the 
entanglement of qubit systems. Here, such bases have natural geometrical meaning. We 
also present 77-kernels of some basic operators acting in A/"-Hilbert space. 

In the A^-Hilbert space of the wave fy-functions one can consider the generalized Schro- 
dinger equation. Formerly, such equation, but in a different form (logarithmic form), was 
proposed by Mandiliara et al. [10]. In the Section 6 we study twocomponent nilpotent 
systems and the eigenvalue problem for qubit system covering various types of couplings 
(Ising, XY etc.). It is interesting that the set of eigenvector for these systems is composed of 
the Bell vectors. There are also considered composed systems of the qubit and boson, the 
qubit and fermion, and the qubit-qubit. As specially interesting we find the qubit-fermion 
system with the supersymmetry acting between two nilpotent parts. 

In the Sec. 7 there is addressed question of the separability of the //-variable dependence 
of the 77-functions and a "duality" of such notion to the entanglement. In the conventional 
function theory, the separability questions are answered using the appropriate Wronskians 
[16]. Here we find that there can be considered special generalizations of Wronski matrices, 
and that their traces and determinants serve as a tool to detect the conditions for separa- 
bility of functions. What is remarkable, one can express relevant invariants known from 
the general theory and entanglement monotones used in recent works on entanglement, in 
terms of Wronskians. This gives additional insight into the structure or flavor of entangle- 
ment. In the description we follow the number of qubits, n = 2, 3, 4. In the present work 
we answer only the question of the pure state entanglement. Finally, conclusions are given 
in Sec. 9. The work is accompanied with several Appendices containing relevant auxiliary 
material and derivations of important facts. 

2. Canonical qubit relations 

As we have mentioned in the introduction we shall follow Wu and Lidar, and by the qubit 
we denote a two-level quantum object with encoded boson-like behavior in many particle 
system, and in the same time, underlying the Pauli exclusion principle. As it was already 
noted in Ref. [12, 13, 14, 15, 4, 10, 8, 5, 9] such a somehow hybrid object is not a fun- 
damental particle, like boson or fermion and inherently carries properties of a composed 
object, but it can be described without any explicit reference to the constituents and the 
way it is composed of. The formalism we develop provides independent description of the 
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qubit itself and gives very convenient characterization of qubit pure state entanglement. 
The nilpotent commuting variables were already used in the context of entanglement by 
Maniliara at al. in [10]. 

In the analogous way as the fermion and fermionic states are understood, one can de- 
fine the qubit system by means of the so called linear commutation relations to have on 
the same footing canonical commutation relations (CCR), canonical anticommutation re- 
lations (CAR) and canonical qubit relations (CQR). To take into account boson-like behav- 
ior of qubit we prefer to use only commutators, and realize the Pauli exclusion principle, 
imposing nilpotency conditions without referring to anticommutators. Such defined object 
is parafermionic, but the form of relations stresses its boson-like properties. This is also 
important from the point of view of supersymmetry. As we shall see later one can define 
supersymmetric system including qubits. Within present scheme, qubits will be even object 
in the graded (supersymmetric) structure, and underly commutation relations as such. 

To describe qubit as an object which is nilpotent but otherwise boson-like, let us ob- 
serve that the commutator of qubit creation operator d + and qubit anihilation operator d 
if nontrivial, cannot have a value in the center of commutator algebra. This would give an 
contradiction with nilpotency. Therefore let us make the following ansatz 

[d,d + ]- = 1 - 2N (2.1) 

On the other hand compatibility condition with the nilpotency of d, d + takes the following 
form 

[d,d + ]+ = 1 + 2Z, (2.2) 
where Z = (d + d — N) is an element from the center of the algebra. Moreover 

[N,d]+=d (2.3) 

[N,d + }+ = d+ (2.4) 
The canonical qubit relations (CQR) consist of the following set of conditions 

[d,d + ]_ = 1 - 2N (2.5) 

[N,d]_=d (2.6) 

[N,d+]_ = -d+ (2.7) 

d 2 = (d+) 2 = (2.8) 

In particular one can take N = d + d. For fermionic operators we have (a + a + ) 2 = 1, and 
here for qubit operators 

{d + d+f = 1 + 2Z, (2.9) 

where for the N = d + d, Z — 0. As we shall show in the following sections, the CQR can 
be naturally realized within the ^-function space by means of //-differential and multipli- 
cation operators. 

Using for the CQR the argument of Wu and Lidar ([4]) one can see that the set of n 
qubits for the large n behaves like a boson. Namely 

1 " 

b = -j=Jldi (2.10) 

v 1=1 
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then 

2 " 

[6,6+]_ = l--y)jV < . (2.12) 

Hence, when the number of sites n is much larger then the effective number of qubits we 

get [&,&+]_ «1. 

3. ^-formalism 

We want to realize the CQR relations in the space of functions of the commuting nilpotent 
variables rj, what suggests the classical limit considered in the Introduction. Because of the 
nilpotency of the fy-variables we can expect some similarities to the superanalysis, but this 
two formalisms are different. The nilpotency is natural for anticommuting (odd) elements 
but for commuting (even) elements it is a result of restrictive condition. In the following 
sections let us collect and introduce some necessary notions. 

(a) Af-numbers 

To introduce necessary functions for describing system of qubits we shall consider 
nilpotent commuting 77-variables defined in Ref. [5] as special elements of an algebra Af. 
We recall, that the algebra N is freely generated by the set of first order nilpotents and the 
unit. Each element v £ Af can be decomposed into the numerical (real or complex) part 
called body and into the nilpotent part called soul 

v = b{v) + s{v) (3.1) 

Any element with a nonzero body has an inverse 

00 

i/- 1 =6(i/)- 1 £(-6(i/)- 1 «(i/)r (3-2) 

m=0 

Analogous formula is known for superalgebra [17, 18]. It is useful to note, that for elements 
with b(v) = 1 it takes the form 

00 

(l + ^i/))" 1 = £(-*(!/))"• (33) 

m=0 

literary generalizing the celebrated formula for real numbers — ^ = ^2(—t) m . 

The /^-variables are first order nilpotents, r\ <G Af , r] 2 = 0. They form subset T> C Af. 
Let r) ^ 77' and t e R, < t < 1 then 

(t V + (1 - t) V ) 2 =2t(t- 1)W (3.4) 

Therefore 77's algebraically dependent with fixed 77 element (i.e. r/r)' — 0) form a star- 
convex set. For algebraically independent elements it is not true unless we admit other, 
higher order nilpotents to enter the set. Further details on the Af algebra can be found in 
[5]. 
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(b) rj-functions 

One can define functions depending on //-variables, where these variables are not fixed 
"generators" but can vary within some set [5]. If we take r)u i = 1, 2, . . . , n to be a set of 
independent nilpotent elements i.e. 

r7f = 0V«, t/i-??2-...t?„^0 (3.5) 

and ff = (?7i, 772, . . . rj n ). Using notation that Iq = and Ik = (ii, ii, ■ ■ ■ ik) is ordered 
multi-index. One can write the expansion of function F(x, ff) € F[ff] of the n ^-variables 
in the general form 

n 

F{x,fJ) = Y J Fi k {x)r 1Ik , (3.6) 

k,I k 

where Fj k (x) E Af and r]i = 1. In the present paper we shall consider only functions 
with Fj k (x) e M, C. This set we shall denote by Fo[x,ff\. Let us observe that there is a 
natural mapping S from the set of symmetric n x n matrices to the ^[r] 1 , . . . , r] n ]. We 
take as components Fj k the principal minors of a matrix with entries indexed by /j. = 
(ii,i 2 ,..., ik), moreover we assume that F@ = 1 e.g. 

B = ( f 11 I 12 ) ~ H(S) - ^V, - 1 + + W + (611622 - 6? 2 )r ? S 2 

V °21 &22 / 

(3.7) 

In the set of 77-functions T n (ff) there exists Z2 gradation related to the decomposition into 
the sets of even and odd functions. The gradation mapping J can be defined as usual by 
relations 

.7(1) = 1, J(%) = -%,and J( V I ») = J(r ] tl )J(r ] ^)...J(r 1 t ») (3.8) 

Another important operation which will be frequently used is duality transformation. It 
resembles the the Hodge * - operator known for exterior forms. The duality operator & we 
shall define as follows 

Qfa/J = Vi n - k , h U I n -k = In, (3.9) 

naturally 0(1) = 771772 • • - T]n and 8 2 = id. Functions such, that Q(F(ff)) = F(ff) we 
shall call selfdual, and these with Q(F(ff)) = —F(ff) antiselfdual. In the sequel we shall 
use notation for this operation or we will denote dualization by *, i.e. *F = Q{F). 

In the next section, in view of further applications we consider some elementary 77- 
functions. 

(i) Elementary ij-functions 

To describe qubit systems we shall use elementary ?7-functions. They can be defined by 
means of a series analogous to the conventional elementary functions. Firstly we introduce 
explicit form of powers of the Tq[x, ff] for n = 1, 2, 3, 4. In the following we shall omit 
the x-variable dependence. Providing explicit formulas we have in mind further concrete 
considerations for n-qubit systems, n = 1, 2, 3, 4. 
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Power function. In the formulas below we sum over strictly ordered configurations 
of indices and appropriate multiplicities of terms are taken into account. 

n=l: 

F( V ) =F + Fit] 

F( V ) m = F™- 1 (F + F 1 r ] ), 
n = 2: 

Ftfrf) = F + Fit] 1 + F 2 r, 2 + F^rfr, 2 = F + Ftf + F^rfrf (3.12) 

Ftfrfr = F^iFo+mF^+F^inFoFu+mim-^F^mm (3.13) 

n = 3: 

Fo + Fit] 1 + F lV 2 + F lV 3 + F 12 r?y + F^rfrf + F 23 rfrf 
+F 123 i 1 1 V 2 r l 3 = F + Ftf + Firfrj* + F ijk r i i rfrf t (3.14) 

F™-\F Q + mFtf) + F™- 2 (mF F 12 
+m(m - l)F 1 F 2 )r ]l m + F™- 2 {mF F 13 
+m(m - 1)^^3)7?!% + F m - 2 (mF F 23 (3.15) 
+m(m - l)F 2 F 3 )mV3 + F™- 3 (m(m - l^Ma 
+F 2 F 13 + F 3 F 12 ) + +m(m - l)(m - 2)F 1 F 2 F 3 
+mF 2 F 123 ) V 1 7 1 2 r 1 3 = F^ l ^ 1 (F + mFtf) 
+F™- 2 {mF F ll + m{m - 1)^)^ + F m - 3 (m(m - 1) 
F FiF jk + m{m - l)(m - 2)F i F j F k + mF 2 F l]k )rfrf 

n = 4: 

F( V \ V 2 , V 3 , r? 4 ) = F + Ftf + Firfrf + F^rfrfrf + F l]kl rf rf rf (3.16) 

F(7 7 1 ,r ? 2 ,r 7 3 ,r ? 4 r - F ( J l ^ 1 (F + mFrf) + F^imF^ + m(m - 1) 
FiFj)ViVj + F ( T~ 3 (™F 2 F ljk 
+m(m - l)F FiF jk + m(m - l)(m - 2) 
■FiFjF k )rfrfr) k + F™-\nF 3 F llkl (3.17) 
+m(m - l)F FiF jk i + m(m - l)F FijF k i 
+m(m - l)(m - 2)(m - 3) • F i F j F k Fi)r] i r] j r] k r] 1 

A cautionary remark: to put the formulas for powers of arbitrary F{ff) in a compact 
form we use the following conventions: when term gets negative power of (F ) m ~ k 
- it vanishes; for the terms with the factor (i*b)° m front, we put (i*b)° = 1 ■ Only 
after such preliminary adjustments we substitute actual values of Fj k . In the present 
work we will be satisfied with the above formulas for n < 4, but generalization to 
higher n is straightforward. 
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Exponent. The r/-exponent we define using conventional expansion, what gives 

e Fw = y EWL = e Fo eS (F m (318) 

Logarithm. Again definition of logarithm function is conventional but it turns out 
that its terms in the expansion have very interesting relations to hyperdeterminants. 
In this context logarithm of functions of commuting nilpotent variables was consid- 
ered in [10]. We shall come back to this in the context of separability. For ry-functions 
with unit body one can define 

ln(l + s(F(ff)) = g(-l)*-i£^M)! (3.19) 
k=i 

where s(F(ff)) is the soul of function F(ff) i.e. here F{ff) = 1 + s(F(rf)). Again 
let us consider explicit formulas for n = 1, 2, 3, 4 using expansions with F = 1. 

n = 1: 

ln(l + s(F(r])) = F\r} (3.20) 
n = 2: 

ln(l + siFirf^ 2 )) = Firf + {Fij - FiFj^rf (3.21) 
n = 3: 

Ml + sfF^ 1 ,!, 2 ,^ 3 )) = Ftf + (Fij - FiFjWr? (3.22) 
+ (F ijk - FiF jk + 2F i F j F k ) V i r ] j r) k 

Using notation f(fj) = ln(l + s(F(r/)) = ln(F(rj)) we can find the following 
identity for components of logarithm for functions belonging to the image of the 
mapping S for n = 3 [19]. Let us note that this identity was derived in [19] without 
any relation to the ry-logarithm. 

A 2 23 = -4/12/13/23 = 46? 2 6? 3 6| 3 (3.23) 
n = 4: 

ln(l + S (F(r/\ r? 2 , ry 3 , V 4 )) = Ftf + (F tj - F^^rf + (F ijk - F t F ok 

+ 2F i F j F k ) V i rfr 1 k + (F m -F ij F kl (3.24) 
- FiF jkl + 2F i F j F kl - 6F i F j F k F,)r i i ifri k v l 

Here as well there exist identities for components of logarithm of function from the 
image of the 5 mapping, namely [19] 

fijkfijlfikl = —4:fiklfijfijfikfjl (3.25) 
^fijklfijfikfil — fikl fijk fij fil "i" fiklfijlfijfik fijl fijk fik fil (3.26) 



Nilpotent quantum mechanics, qubits, and flavors of entanglement 



11 



;onometric functions. This family of ^-functions is defined by formal series 
ogous to the conventional one, namely let F = Fi k V Ik > Fi k € ^ 

C0S (F(r?)) = £(-l) fe -— (3.27) 



Trig. 

analogous to 



(2k) 

p2k+ 



sm (F(r i ))=J2(-l) k WTT y^ 0.28) 



fc=0 



Obviously cos 2 (F) + sin 2 (F) = 1 and we have 

cos(F) = cos(F )cos(s(F)) - sin(F )sin(s(F)) (3.29) 

sin(F) = sin(F )cos(s(F)) + cos(F )sin(s(F)) (3.30) 

for a function of one n variable we obtain that 

cos(F(r])) \ = ( cos{F ) -sin{F ) \ f 1 
sin(F(r))) J \ srn(F a ) cos(F Q ) J \ F lV 



(3.31) 

For further considerations we will need sin(J2rf) and cosi^n 1 ) in explicit form for 

n = 2,3,4. 

n = 2: 

costf+rf) = 1 77 1 ?7 2 (3.32) 
sin^ 1 + n 2 ) = n 1 +n 2 (3.33) 

Let us note that cos(n 4 + rj 2 ) is antiselfdual and sin^ 1 + rj 2 ) is selfdual. 

n = 3: 

cos(rj 4 + rj 2 + -q 3 ) = 1 - 7?V - 77V - f?V (3-34) 
sm(r) 1 + if + rf) = n 1 + rj 2 + rf - n l rj 2 rf (3.35) 

Here 6(cos(r? 1 + rj 2 + r/ 3 )) = -sin^ 1 + nf + rj 3 ). 



d = 4: 

cos (77 1 + r/ 2 + 77 s + rj 4 ) = I - rfn 2 - rj l n 3 - rfrj 4 - tfrj 3 (3.36) 

2 4 3 4, 1234 

— -q -q — -q rj + rj rj rj rj 

sinirj 1 + rj 2 + ij 3 + rj A ) = rj 1 + n 2 + rj 3 + rj 4 - rj 1 rj 2 l rj 3 - rj 1 rj 2 V (3.37) 

- rfrf-q 4 - rfrf-q 4 

In this case cos(rj x + rj 2 + n 3 + n 4 ) is selfdual and sm(rj 1 + n 2 + rj 3 + rj 4 ) is antiselfdual. 

(c) Distinguished rj-functions 

It turns out that //-functions be related to combinatorics and graph theory. Here we 
want to describe a family of rj functions yielding the normal and ordered Bell numbers and 
Stirling numbers of the second kind. 
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Functions: T n and £ n . Let us define the following fy-functions 

n 

% l = T n (r h ,r l2 ,...,r] n )=J2v i (3-38) 

i=l 

and 

n 

£n = £n{Vl,V2,---,Vn)= ^ (339) 

fe=l,J fc 

As we shall see they play important role in describing combinatorial characteristic 
in terms of 77 variables. Directly from definitions we have the following relations 

£ n = e Tn - 1 = s(e T ") (3.40) 

Expanding k th power of £ n we find that numerical coefficients in respective terms of ex- 
pansion just count the number of functions from m-element set into an fc-element set, 
giving the following formula [20] 

i £ n) k = £ 5>!S(m,fc)»7 7m , (3.41) 

k<m<n I m 

where S(m, k) are Stirling numbers of the second kind. Because the following sum gives 
the Bell number B n (so called normal Bell number) 

B n = S ( n > m )' 71 ^ (3.42) 
from the Eq. (3.41) we get that 

e '" = EE B ^ 4 ' (3 - 43) 

k=0 I k 

It is easy to prove another important relation which involves ordered Bell numbers C n , 
namely 

k=a i k 

For convenience of the reader let us recall the values of Bell numbers for n = 0,1, ... ,5. 
Namely, B n : 1, 1, 2, 5, 15, 52 and C„ : 1, 1, 3, 13, 75, 541, respectively. 

(d) Symmetric functions 

Conventional symmetric polynomials are ubiquitous in many areas of mathematics and 
mathematical physics. Specially their relation to the representation theory and theory of 
invariants makes them to appear in many physical applications (to name a few: Boson- 
Fermion correspondence,integrable systems, super/string theory, Chern-Simons theory, link 
invariants and 3-manifolds invariants versus moduli spaces of Riemann surfaces). Here we 
want to adapt conventional theory to the case of nilpotent commuting 77-variables, having 
in view the applications to multiqubit systems. 

We shall call a F{ff) the symmetric ^-function (symmetric ^-polynomial) if 



F{Va(l),Va(2),---,Va(n)) = F(Vl,m,---,Vn), O G S n (3.45) 
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Then we can consider the elementary symmetric ^-polynomials e k = e k (if) 

e = 1 (3.46) 

n 

ei = ^m = T n (3.47) 



e2 



= £wj (3-48) 



(3.49) 

e„ = ?7i772 ■■■Vn = £n (3-50) 
From above definition and formulas for powers of F(ff) we have relations 

e k = k\ek, or eiefc_i = fce4 (3.51) 

moreover O(ej) = e„_j. Hence, for the even n there exists selfdual polynomial eq. The 
set of symmetric ?7-polynomials is naturally Z„ -graded, by the degree of polynomial. 
An arbitrary symmetric ry-polynomial s k has expansion Sk(ifj = / s h r lii- An s k is 
homogenous of degree k if Sh(Xrf) = X k Sk{fJ), where A e M. As in the conventional 
case we can introduce the analog of Euler operator counting the degree of the homogenous 
polynomial. Namely, 

ff-V = Y,Vidi, (3.52) 

i 

where di denotes derivative with respect to the rji variable ([5]) (cf. also the next section). 
In the context of nilpotent quantum mechanics we shall call this mapping - the qubit num- 
ber operator. 

Let us formulate the fundamental theorem of symmetric ^-polynomials: the F(ff) is sym- 
metric polynomial iff it can expanded in elementary symmetric polynomials 

n n 1 

F(rJ) = F(ei, e 2 , . . . , e„) = £ F ]k] e k = ^ T^|fc|^ (3.53) 

k=0 k=0 ' 

In analogy to conventional theory we can define the complete symmetric ^-polynomials 

hk(ff)= J2 Vf^---Vi n , (3.54) 

dl+d 2 H Yd n =k 

where di = 0, 1. Using real parameter t they are generated by the function 

H n (t) = ]T h k (ff)t k = Yl ^V? ■ ■ ■ V d n H d ^+-+ d " = 1 

fe>0 d u d 2 ,...,dn Hi=nU t7 li) 

(3.55) 

Analogously, generating function of elementary symmetric ry-polynomial ek(ff) is given as 

n n n 

E n {t) = ]T e k(v)t k = J2 T\ e *W* = exp ^ + ^ + •■■»?») = IT ^ (3 - 56 1 

k=0 k=0 ' k=l 

Therefore 

E n (-t)H n (t) = 1. (3.57) 
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Hence, for m > 1 



^2(-iye k h m ^ = 



(3.58) 



fc=0 



As in conventional case one can write above relations as vanishing determinants. 
Examples: 

• m=l. eo/ii — ei/io = 

• m=2. 

/i 2 = - e 2 = 

and hence 2e 2 = e\ 
For arbitrary /c, the identities for the e k now read as 



ei e 2 
1 ei 



ei 


e 2 • 


■ e fe _i e fe 


1 


ei . 


e k -i 





1 . 


e k -2 







1 ei 



(3.59) 



Finally let us observe, that the notion of antisymmetric ^-polynomial a„(?7 (T ( 1 ) , r? CT (2),..., ^(n) 
(— l)l cr la n (?7i, r/2, • • • , T) n ) is almost trivial. Namely, the Vandermonde determinant A(rji, ?y 2 , 
n, <i<j<nijli ~ for nilpotent commuting variables can be different form zero only for 
n = 2. 

(e) r\-calculus 

^-derivative. The ^-derivative in form we need in the present approach was pre- 



sented in Ref.[5]. To fix the notation let us recall its main properties. Namely, let 

dj = £ (3.60) 

then 



dirf = Sf, d t l = 0, didj = djdi 



(3.61) 



Instead of the Leibniz rule for F(ff), G(fJ) £ F[ff\ we have the following relation 
(para-Leibniz roule [6]) 

d l (F ■ G) = d t F ■ G + G ■ 8,F - 2 Vi d i Fd i G (3.62) 

For further details cf. Ref. [5, 9]. 

^-integration. The ^-integral is defined by the following contractions on basics vari- 

(3.63) 



ables 



J V l dv 3 =5), J d m 







and by linearity is extended to the Fo [ff\ . Despite the fact that definition is exactly the 
same as for the Berezin integral, the properties of the multiplication in the algebra 
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N make that in this formalism the integration by part formula has different form. 
Namely, 

J Fdr^J ^ Gdr^j = X - (^J (dF) ■ Gdr\ + J F ■ (dG)dr^j . (3.64) 

Detailed properties of this integral are given in Ref.[5]. Below we recall Gaussian 
integral (firstly considered by Palumbo [12]) which we shall use to define scalar 
product in F [ff\. Then we introduce the 77-Fourier transform. Finally we show the 
form of generating //-functions for Stirling numbers of second kind and the 77-integral 
form of the triangular recurrence relation for S(n, k). Moreover we give //-integral 
expression for both types of Bell numbers. Because in the definition of the 77-Hermite 
polynomials we make use of the notion of //-derivative, therefore we introduce it in 
this paragraph. 

Gaussian //-integral. Let Bbeanxn symmetric matrix then 

^'drjdrj' = per(B), (3.65) 



where per(B) is the permanent of the square matrix B (cf. Appendix A) and dff = 

drfn = dt]idt]2 . . . drin 

/y-Fourier transformation. Let us define 77-exponent as follows 

e <lv> = Y^^<lri> k , <£,r?>=f%, (f) 2 -fe) 2 =0. (3.66) 

k 

A Fourier //-transform of a function g(ff) e F[ff\ we shall call the function g(t;) G 
such that 

(^)(D - 9(0 = J e<^>g(ff)drf. (3.67) 
Then the inverse Fourier /y-transform is given in the following form 

T-\f(rj) = J (3.68) 

Such defined Fourier ?7-transform has properties similar to the Fourier transform 
considered in superanalysis. For example 

%(^)K) = ^fe)K) (3-69) 

Hliv Va)M) = e<^>(ff)(0 (3.70) 

+ &)) = ne <c - ,J> fm 0.71) 

There is also interesting property of Fourier //-transform obtained when we want to 
transform the //-coordinates. Let A be invertible n x n matrix with entries from Af 
then realizing permutation and scaling transformation 

HArf)(i) = Tf(A- 1T i) ■ per(A), (3.72) 
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In opposite to the case of the Grassmannian Fourier transform, we do not obtain 
here the determinant but permanent. This is characteristic property of 77-calculus. In 
particular, for n = 2, 4 we have that 

T(Am) = Tf{A- lT i) ■ Hf{A), (3.73) 

(cf. Appendix A.). 

77-integral form of the Stirling numbers. Using directly the expansion of £ n we 
can represent the Stirling numbers of the second kind by the following //-integral 

S(n, k) = ^J (£n) k drf n (3.74) 

It is easy to see that triangular recurrence relation for S(n, k) 

S(n, k) = k S(n - 1, k) + S(n - 1, k - 1) (3.75) 
can be expressed by means the following compact ^-integral relation 

J(£ n ) k drf n ^k JiS^-'e^drfn-i (3.76) 

77-integral form of the Bell numbers. Such form of the normal Bell numbers is di- 
rect consequence of the expansion (3.43) and was firstly given in [21] 

B n = J e £ -drf n (3.77) 

Here we obtain new formula for ordered Bell numbers. Using (3.44) and relation 

2 - e r " = 1 - £ n (3.78) 

we get 

C n = J (2-e r ")- 1 dr7„ (3.79) 
The last formula resembles relevant conventional representation of the Bell numbers. 

77-Hermite polynomials. Let <C fj, ff »= \e\(ff). By the 7 fe -th 77-Hermite polyno- 
mial we shall understand function of the form 

H Ik = ( m ,m, ■■■,Vn) = (-l) k e«^»d Ik e-«^, (3.80) 

where the multi-index Ik is fixed. Now, the fc-th degree 77-Hermite polynomial is 

H\k\ = X! Hlk > (3.81) 
ik 

where the length k of the multi-indices Ik is fixed. are symmetric polynomials 
e.g. for n = 2: i?| | = 1, = e\, H\ 2 \ = — e — e\\ for n = 3: H\ a \ = 1, 
H\i\ = 2ei + 6e 3 , H\ 2 \ = -3e - 3e 2 , H\ 3 \ = 0. 
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4. Qubit systems in Af-Hilbert space 

To describe qubit and many-qubit states let us use the structure of bimodule over N alge- 
bra. Let H be such Af-module. We equip it with the TV-scalar product i.e. .AAsesquilinear 
mapping 

<.,.>:HxH^N (4.1) 

such that for F,G eU 

<vF,G> = <F,vG>, v&N (4.2) 
< F,G> = VG€H=>F = (4.3) 
b(<F,G>)* = b(<G,F>) (4.4) 
&(< F, F >) > 0, VFeH (4.5) 

(4.6) 

We shall call such module with TV-scalar product the A^-Hilbert space. Vectors F with 
&(< F,F >) 7^ we shall call physical. Particular realization of such TV-module is 
given by the function space F[ff n ]. In this module we introduce TV- valued weakly non- 
degenerated scalar product in the following form 

< F, G > N = J F*(rf)G(ff)e<^^> drf* dff, = J F* (rf)G(rf)d{i(rf , rf) (4.7) 

where 

n 

^=EE f ?y** (4 - 8) 

k=0 I k 

and ★ denotes complex conjugation. For FoiVn] obviously Fj k are complex numbers. The 
first order nilpotents rf* are algebraically independent from rf. In components we have 

<f, g> n =y,Y. f i g i« < 4 - 9 > 

k=0 i k 

To see how 1 -qubit algebra is realized in this formalism let us take T[rf\ (set of ^-functions 
of one variable). In particular, ^-scalar product of F(rj) and G{rf) functions takes simple 
form 

< F, G >aa— F* Go + F*d (4.10) 
and in this space realization of the qubit algebra is given in the following form 

d + =Tj-, d^d/drj (4.11) 

i.e. operators d + and d are conjugated with respect to the scalar product given by (4.9) 
moreover 

a 3 = l-2rfd v (4.12) 

is self-conjugated and [d, d + ] _ = a 3 . Hence, we obtain natural realization of the canonical 
qubit relations Eq.(2.5-2.8), with N = rfd v and [TV, d]_ = d, [TV, = -d + . 

Despite the scalar product we can define in the n-qubit space another weakly non- 
degenerate form, symmetry of which depends on the parity of n. Namely, the gradation 
mapping J allows to consider the natural orthogonal projections 

7T ± F(rf) = ^(F(rf)±J(F(ff))) (4.13) 
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on the even F + and odd F_ part of the function F. Using the mapping J we can define the 
following linear weakly non-degenerate form 

Lo n {F, G)= f J(F(rf)yG(fJ)d m ...d Vn = ]T(-l) fc f/ fc G/„_ fc , where I k Ul n - k = I n 

J k,I k 

(4.14) 

The n is symmetric or antisymmetric, depending on the parity of n 

uj n (F, G) = (-l) n w„(G, F) (4.15) 

Let us note here that it is the 77-version of the form obtained in the tensor product of C 2 ® 
• • • ® C 2 from the antisymmetric form e in the C 2 . Explicitly, we have 

n=l: 

G) = F G 1 -F 1 G (4.16) 

n=2: 

u) 2 (F, G) = F G 12 + F 12 G - F X G 2 - F 2 G X (4.17) 
In particular for basis {1,7/1, V2, I1I2} we get 

^2(1, 771T/2) = 1, W2(r?i,r? 2 ) = 1 (4.18) 

and in this case we have exactly the form considered by Wallach [22] 

n=3: 

u>3(F, G) = F0G123+F123G0+F23G1—F1G23+F13G2—F2G13+F12G3—F3G12 

(4.19) 

Counterpart of the above symplectic form was discussed by Meyer and Wallach 
[23], where identification of components (xi,yj) i,j — 1,2,3,4 used there and 
components of the ^-function is the following: xi — F , x 2 — F23, x 3 — F13, 
Xi = F 12 , yi = F 12 3, V2 = F lt y 3 = F 2 , y 4 = F 3 . 

In the space of ry-functions we can introduce natural family of projectors related to the 
decomposition of the F(ij) into the part depending on the fixed and independent of 
it, F(jf) = F(7]!, 772, • • • , fjk, ■ ■ ■ ,f]n) + VkFitj!,^, r] n ), where hat indicates 
skipped variable. Natural realization of it is given by 77-derivative operator, due to the iden- 
tity 

dkVk + Vkd k = 1 (4.20) 

we can introduce projectors 

7Tfe|o = d^k ■ (4.21) 
7Tfc|i = Vkd k - (4.22) 

For a fixed k they are orthogonal and for different indices k they commute 

^k\i^k\j = SijTT k \j, 7T k \0 © 7Tfe|l = id (4.23) 

Kk\m\j = n\j^k\i, k ^ I (4.24) 
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For example the decomposition of the F(rji , 772) has the following form 

F{m,m) = F + F 2 f]2 + Vi(Fi+F 12m )=F( m ) + mF(V2) (4.25) 

= ^i\oF(Vi, V2) +7Tl|l^(7?l, %) 

or 

F(m,m) = F + F iril + r, 2 {F 2 + F 12 rn) = F (771) + 772^(771) (4-26) 

= ^lo-FOft, + 7r 2 |iF(r7i, 772) 

Moreover, we get the full decomposition of the ^(771 , 772) using composition of projectors 

^0 = ^i\i^2\oF(r]i, 772), (4.27) 
Fim = ^111^210^(^1, m), (4-28) 
F 12 mV2 = ^i\i^2\iF(m, V2) (4.29) 

In the naturay way this formula generalizes to the n variables. 

For n = 1 we can realize the antisymmetric form u>i in alternative way by introducing 
the wedge product of 77-functions as 

F(r?) A G(r?) = (F Gi - F\Go)rj, (4.30) 

then the anti-symmetric form obtained in this way we shall denote T>i 

T>x{F, G) = J F(rj) A G{r,)d v (4.31) 

For one 77 variable T) 1 (F, G) = u>\. It is possible to generalize the T>i to the d — 2 case, 
in the following sense. The space of 77-functions F(r]i, rj 2 ) is 2 2 -dimensional, so admits 
antisymmetric form. Because, a general d = 2 77-function can be decomposed using one of 
the factorizations (4.25) or (4.26) then we can define entity C 2 , which is counterpart of the 
concurrence of the 2-qubit states, as 

C 2 = P 1 (F(77 1 ),F(t /1 )) = P 1 ( J F(772),F(77 2 )) - F F 12 - F^ 2 (4.32) 

The anti-symmetric form V 2 we define in as follows 

V 2 = V 1 (F( m ),G( Vl ))-V 1 (G( m ),F( m )) (4.33) 
= ViiFfo), 6(772)) - Pi(G(r7 2 ), F(t? 2 )) 

(a) Natural bases: monomial, trigonometric 

Having in mind further application to the n = 2, 3, 4 qubit pure states entanglement let 
us discuss particular bases in the !F[ff n ]. 

(i) monomial basis 

The simplest basis in the space of 77-functions is the monomial one consisting of 
{vi k }k=o- I* correspond to the so called computational basis in conventional notation, 
widely used in the literature. The binary notation used there for describing the elements 
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of the tensor product of C 2 is simply related to the 77-notation. Namely, translation of the 
multi-index for system of n qubits with the "binary" entries 0, 1 to the multi-index Ik used 
in the //-function expansion is obtained by putting ordinal numbers equal the position of l's 
appearing in the binary multi-index e.g. (0, 0, 0, 0) 1— > 0, (1, 0, 0, 0) 1— » 1, (0, 1, 0, 0) 2, 
.... (0, 1, 0, 1) 1 ^ (2,4), (1, 1,1, 1) 1 ► (1,2,3,4) (cf. also Appendix B). Using such 
notation we obtain for example the following expression for the hyperdeterminant of the 
hyper-matrix B = (bijk), i, j, k = 0, 1 

Det(F) = Det(B) = (F*F? 23 + F*F? 2 + F*F? 3 + F^F^) + A(F F 23 F 13 F 12 

+ FiF 2 F 3 Fi 23 ) — 2(F F 3 Fi 2 Fi 23 + F a F 2 F 13 F 123 

+ FqFiF 23 Fi 23 + F 2 F 3 F\ 3 Fi2 + FiF 3 F 23 F\ 2 

+ F 1 F 2 F 23 F 13 ) (4.34) 



(ii) trigonometric basis 

More interesting from the point of view of the entanglement, there are trigonometric 
bases. For the n — 1 such basis is identical with the monomial one {1, 77}, but for higher n 
we obtain nontrivial functions. 

n=2: We take trigonometric function (3.32) with arguments 771 ± rj 2 . These functions 
are orthogonal and when normalized with respect to our A^-scalar product they take 
the following form 



hi = 


1 r 1 

— = COS(77 

V2 


+ v 2 ) = 






(4.35) 


h 2 = 


1 f 1 
— = cos(77 

V2 


-v 2 ) = 


> + *V> 


= ?pGHZ+ 


(4.36) 


h 3 = 


1 ■ ( 1 
—= sinvn 

V2 


+ v 2 ) = 




= ipw+, 


(4.37) 


hi = 


1 ■ t 1 
—= svniri 

V2 


-v 2 ) = 




= tpw-, 


(4.38) 



One can recognize here the 77-realization of the "magic" basis for 2-qubit system. The 
cos-states are GHZ-type and sin-states are W-type (Bell states up to the particular 
phases). According to our definition, above states are physical with respect to the 
introduced before 77-scalar product with "Gaussian measure". 

n > 3: Above basis can be generalized to higher n, it contains 2™ elements 

f 1 1 

: C0S(771 ± 77 2 ± 77 3 • • • ± Vn), 7= = Sm(77l ± 772 ± 773 • • • ± T] n ] 



{ a/2"- 1 V2™- 1 

(4.39) 

Functions entering above basis are normalized. 

(b) rj-kernels 

We can realize the action of an operator in A^-Hilbert spaces of 77-functions using "ker- 
nels" with respect to the 77-integral. 



J A(rf,rj 



(AF)(rf) = / A(ff,ff)F(ff)drf (4.40) 
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Taking expansion of the kernel in the following general form 



A(ff,ff) = J2A IklJl 



(4.41) 



we obtain the following explicit action on components of the function F 

(AF) Ik =J2 A ^Fj n -n (4-42) 
i,Ji 

where J„_; is complementary strictly ordered multi-index i.e. J/ U J n -i = Jn- In the 
fundamental d = 1 case we have 

AF{rj) = (.Aoio-Fi + A oll F ) + (AnoFi + A M1 F )ri. (4.43) 

and one can easily obtain realizations of some operators important in Hilbert 77-space 



id 


= A id (r], r)') 


= 11 + 1]' 


= S( V - il') 


= stf - v) 


(4.44) 


V 


= A d+ (r],r]') 


= r\rj 






(4.45) 


d v - 


= A d (r], rj) -- 


= 1 






(4.46) 


d v V 


= A 7 r (r],r)') 


= r{ 






(4.47) 


r]d n - 


= A Wl (r],r]') 


= V 






(4.48) 


a 3 - 


= A a3 (r), 77') 


= r{ - 7], 


cr 3 = 1 - 


2r/d v 


(4.49) 


ay 


= A ai (r],r]') 


= e^'i = 


cos(rj — rj 1 ) 


= ch(T] + T]'), 


0-1 =d v +ri (4.50) 


e- 


= A e (r,, V ')-- 


= e-* 1 ' 1 ! = 


cos{rj + rj') 


= ch(7]-T]'), 


(4.51) 


e 


= —ioi = d v 


- V, 









were symbols Oi, i = 1, 2, 3 are used because above 77 operators play the role of c-matrices 
in Tj-realization of the su(2) algebra and ty/,, k = 0, 1 is the projection on the first and 
respectively on the second term of the F(r/) expansion. The Hadamard operator (gate), 
defined by relations H 1 = -^=(1 + 77) and_ffr7= -^(1 — 77) has the following 77-differential 
realization 

1 



H= -^=(l + 77 + <9 t? - 277<9„) 



(4.52) 



and its ^-integral kernel is of the form 



1 



H- = A H (r], rj 1 ) = -j=(l-ri+rj' +rii] 1 ) = ch(r]' +r])+sh(r]' -r]) = cos(t]' -r])+sin(r]' -r]) 

(4.53) 



5. Nilpotent quantum mechanics and 77-Schrodinger equation 

As it is well known from supersymmetric theories, we can realize classically fermions us- 
ing anticommuting variables and formalism of supermechanics or pseudomechanics. There 
we have prequantum description of e.g. spin systems, using anticommuting ^-variables and 
using explicitly supersymmetry or not. The superphase space language, graded Hamiltoni- 
ans, graded Poisson brackets, graded Heisenberg group, graded special functions etc. ap- 
pear very useful in classification and analysis of properties of such systems. Then applying 
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cannonical quantization we obtain quantum description of such systems. With this respect 
there are two approaches. One procedure yields after quantization the conventional Hilbert 
space formalism and conventional quantum mechanics. Here symbolically the Grassmann 
algebra of ^-variables is traded for Clifford algebra of er-matrices (or 7-matrices in rel- 
ativistic case). It is so called Casalbuoni/Berezin=Marinov quantization ([24, 25]). But 
there is also another version of quantization procedure, which preserves the ^-variables. 
This one is specially efective when we ask about representations of supersymmetry, look 
for multiplets of states and want to use the Feynman path integral for fermions. It can be 
symbolically named the super-Schrodinger quantization, because we use there superwave- 
functions i[)(x, 9) and generalized Schrodinger equation, involving superderivatives. 

In description of qubits we want to develop analogous approach to the super-Schro- 
dinger one. Some elements of it are already known in literature. On the one hand, in the 
papers by Mandiliara at al. [10, 11] there was already used equation which can be named 
?7-Schr6dinger equation in the logarithmic form. It was written for a restricted set of func- 
tions of nilpotent commuting variables - nilpotentials, and used to address, via system con- 
trol methods, some questions of the entanglement. On the other hand there exists classical 
theory based on nilpotent commuting variables introduced in [5] which provides config- 
uration and phase space description of nilpotent systems. It is called nilpotent classical 
mechanics. Moreover another essential aspect of such theory i.e. path integral formalism, 
was discussed some time ago by Palumbo at al. [12, 13, 14]. Let us note that in all this 
approaches, except the nilpotent classical mechanics there was neglected the fact, that the 
derivative with respect to the nilpotent commuting variables do not satisfy the Leibniz 
rule, what makes the whole construction nontrivial. So, it is natural to consider the nilpo- 
tent quantum mechanics as formalism which is related by a "^-canonical quantization" to 
the classical nilpotent mechanics. Because known by now ?7-Poisson brackets do not sat- 
isfy the Jacobi identity, the term "^-canonical quantization" leaves some open questions, 
but the formalism of nilpotent quantum mechanics itself is consistent and effective. We 
shall use here the restricted ?7-Schr6dinger quantization in the following sense. To quantize 
classical nilpotent system, we take a classical observable in the normal ordered form i.e. 
momentum variables are to the right of the coordinate variables and realize position and 
momentum as operators 

Q 

Vk — >f) = Vk-, P — >Pk = ^— (5-1) 

in the 7V-Hilbert space of 77-functions depending on t)k, k = 1, 2, . . . , n. Let ip € T[x, ff[ 

ih^ip(x, rj, t) = Hip(x, rf, t), (5.2) 

where H is quantized Hamiltonian H(x,p x ,r],p ri , t) of the system. For the two level sys- 
tems it is typical to consider explicit time dependence of the Hamiltonian. For example in 
the n — 1 case the Hamilton function is singular in its nilpotent part in the sense that it 
contains terms linear in p v i.e. H = -^p x + b(t)p v + c(t)i]p ri + V(x, ij, t). After quantiza- 
tion we can write this Hamiltonian in the convenient form H = ^p x + V(x) + B{t) ■ a, 
where nilpotent part can be written as 

d d 

Hnil P = (B x (t) + iB y (t)) V + (B x (t) - iB y (t))— - 2^(t)r7— + B z . (5.3) 
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In the present paper we restrict ourselves and analyze the properties of the nilpotent part 
alone, putting aside the question of the simultaneous x coordinate dependence. We will 
assume as well the global factorization of time dependence of the r\ - wavefunction il>(ff, t) 
and study the stationary ?7-Schr6dinger equation for nilpotent quantum system. 

H^ff) = \^(rf) (5.4) 

The structure of eigenstates for multiqubit systems turns out to bo nontrivial, when one 
addresses the question of entanglement. 

Let us compare considered above the 77-Schrodinger equation with the one studied by 
Mandiliara et all [10]. In the latter one, the authors restrict to the case when 77 - wave 
function has invertible values in the algebra TV (we describe situation using formalism 
developed in the present work) i.e. -0(r/) = tp + ij)^ + . . . and ip ^ therefore one can 
take function tp(ff) = ^i'i'n) an d there exists its logarithm f(ff) = Inip. Now, because 

ij t f(rf) = ij t lni>(ff) = z^^H^rf) (5.5) 
and ip(ff) — e^v) we g e t 

ij t f{ff) = e- f ^He f ^\ (5.6) 

what is the form of the equation employed in [10]. But let us note once again, there many 
other states that have r\ - wave functions with non-invertible values (like Werner -like states 
) and for them such equation is not valid, but fy-Schrodinger equation (5.2) can be used 
without problem. 

6. Two component nilpotent systems 

Let us consider well known generic Hamiltonian for the two, two level quantum systems. 
In terms of a - matrices it is given by 

H ClC2 c 3 = cicr x ® a x + c 2 a v ® cr v + c 3 a z a z , (6.1) 

where c,, i = 1, 2, 3 are numerical parameters. In the ?7-Schr6dinger representation it can 
be written in the following form 

H = ( Cl -c 2 )(d + ®d + + d<g>d) + (c 1 +c 2 )(d + ®d + d®d + ) (6.2) 
- 2c z {d + d® 1 + 1 ® d + d) + 4c 3 (d+d«) d + d + i), 

where as before d + = r\- and d = d, q . There are several special choices of the values of the 
parameters discussed in the literature: 

1. Ising coupling: c = c\ ^ 0, c 2 = c 3 = 

His = c(d + ® d + + d <g) d) + c(d + (g>d + d(g>d + ) (6.3) 

It is interesting to note, that classical Hamiltonian for this system is composed of r/- 
harmonic oscillator and additional angular momentum like term i.e. H = c(pip 2 + 

mm + V1V2 + mvi) 
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2. XY coupling: c = c\ = c 2 ^ 0, C3 = 

Hxy = Md + ®d + d®d + ) (6.4) 
Here the ^-harmonic oscillator part is not present 

3. Heisenberg (spin system) type: c = c\ = c 2 , c 3 7^ 

H He = 2c(d + ®d+d®d + )-2c 5 {d + d + ®l+d®d + d)+Ac 3 {d + d®d + d+^) (6.5) 

Here there it is a subcase with c = 0. 

4. Q-invariant type: c = c\ = —c 2 

H Q = 2c(d + ®d + +d®d)-2c 3 (d+d + ®l + d®d+d)+4:C 3 (d + d®d + d+j) (6.6) 

Before solving the 77-eigenfunction problem, let us observe that for the H ClC2C3 family of 
Hamiltonians there exist conserved charges. One is related to the duality transformation O. 
For the two-qubit system the is implemented by operator 

e (2) = d®d + d+ ®d + d®d+ + d + ®d + . (6.7) 

We have that 

[6< 2 \ i? ClC2C3 ]- - 0, VC1,C2,C3 (6.8) 

Let us define two charges 

Q = d®d+ (6.9) 

Q + =d + ®d (6.10) 

The Q and Q + are conserved only for the H ClC2C3 with c\ + c 2 = 0. 

Let us now study in the ry-Schrodinger formalism, the structure of the eigenstates of 
this composite system. The H ClC2C3 Hamiltonian has explicitly the following realization in 

H ClC2C3 = ((ci + c 2 )i]2 - 2c 3 ij 1 )d 1 + ((ci + c 2 )m - 2c 3 i] 2 )d 2 (6.11) 
+ (ci - c 2 )??i??2 + ((ci - c 2 ) + 4c 3 77i7?2)9i5 2 + c 3 

For the eigenproblem 

H ci c 2 c 3 4>(vi, m) = W(vi, m) (6- 12) 

the set of solutions is shown in the Table 2. For c\ — c 2 = there is degeneration with 
Ai = A 2 and analogously A3 = A 4 for c\ + c 2 = 0. Because the duality transformation 
commutes with the Hamiltonian the eigenvectors have to be selfdual or antiselfdual i.e 
0(ipi) = tpi, 6(?/> 2 ) = — tpi, 6(^3) = —"03 an d 6(^4) = ^4- A nontrivial action of 
duality transformation is in eigenspaces with degenerate A. 

Besides the considered above qubit-qubit system, we shall discuss the system with 
supersymmetry which is composed of qubit and fermion (Q-F) in analogy to the standard 
case of boson-fermion (B-F) system. In addition we consider below the boson-qubit (B-Q) 
and fermion-fermion (F-F) systems for which at least one part of the composed system is 
nilpotent. In these cases one can add to the conventional SUSY-like Hamiltonian the Forster 
term, which gives exchange interaction and is characteristic for the two level systems [26], 
but breaks the SUSY. 
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Table 2. Eigensy stem for n = 2 (tpo, i>i, i>2, f>i2 are arbitrary). 



Ai 


= Ci - C 2 + C 3 


*1 
*1 


= 75(1 + ^) 

= IpO + V>12?)1»?2 


ci - c 2 j= 
ci — c 2 = 


A 2 


= ~(C1 - C 2 ) + C 3 


* 2 
* 2 


= 75( 1 ~ ? ? 1 ^) 

= IpO + 1pl2VlV2 


ci - c 2 / 
ci — c 2 = 


A 3 


= -(Cl + C 2 ) - C 3 


*3 
*3 


= -72^1 

= IplTjl + tp2T]2 


ci + c 2 / 
ci + c 2 = 


A 4 


= Ci + C 2 — C3 


*4 
*4 


= -72^+^) 
= IplTjl + 1p2V2 


ci + c 2 / 
ci + c 2 = 



(a) boson-fermion 

To fix the notation let us recall the well known SUSY system of quantum bosonic and 
fermionic oscillators. Hamiltonian of such a composed system has simple form 

ff< ) = ui (1 ® /+/ + ® 1) = ^1 ® (/+/ - 1) + (&+& + ^) ® l) (6-13) 

It is a special case of more general Jaynnes-Cummings Hamiltonian which is supersym- 
metric only for special values of parameters. Namely, for u — dj an d K = 

Hj C = lo (b+b - i) ® 1 + c^o 1 ® (/+/ - 1) + «ff (F) , (6.14) 

where i?( F ) is the Forster Hamiltonian [26] of the form 

= b+ ® / + 6 ® /+ (6.15) 

To have consistent grading the H^ F > and k parameter are odd entities. Let us write ex- 
plicitly SUSY algebra for the Jaynnes-Cummings system. The supercharges are defined 
as 

Q B F = iVub®f + , Qg F = —iy/uj b + ® /, (6.16) 
and hamiltonian has the form 

H ( ^=QI f + Qbf, H BF = Hf F + nH { ^ (6.17) 
Graded commutation relations now read as 

[Q + BF , Qbf}+ = Hf F , (6.18) 

[Hbf, Q + bf ]- = kH™, [Hbf, Qbf}- = kH™ (6.19) 

(b) qubit-fermion 

For this system composed of qubit and fermion we have the odd supercharges as in the 
conventional boson-fermion case 

Q QF = iy/U d ® /+, Q+ = -iy/u d + ® / (6.20) 
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Graded commutation relations now read as 

[Qqf, Qqf}+ = w (d+d ® 1 + (1 - 2N d ) ® /+/) = Hqp (6.21) 

H ( qf = Qqf + Qqf, H qf = H§1 + kH$ (6.22) 

[Hqf, Qqf}- = kH™, [Hqf, Q% f ]- = nH Q % (6.23) 

This system is very interesting because both of its parts are two level systems and we 
still have nontrivial supersymmetry transformations with graded algebra of charges (for 

H Q F F = 0)- 

The Hamiltonian Hq f in Schrodinger representation takes the form 

H ( Q % = cj( v d ri + 9d e - 2 v 6d v d 9 ) (6.24) 

Solving generalized stationary 77-Schrodinger equation Hg F F(r], 6) — XF(rj, 6) we get 
that there is eigenspace related to the zero energy, and another one to A = u. So the 
vacuum is degenerated and invariant under supersymmetry transformations. The subspace 
with nonzero energy is degenerated, as it should be in supersymmetric system. The non- 
unique ground state is peculiar. So, we have cf>o = 1 and ip — r/9 even and odd respec- 
tively, ground states (in the sense of the Grassmannian parity) and <f> u = t], = 9 even 
and odd excited states. 

H§l<f> u - u>4> u , H$ll> u =unl> u , ff$.0o = O, H Q a) F yj = 0; (6.25) 
= *l> u , - 0, Q(j) - 0, Qyj = (6.26) 

Q + 4>u = 0, Q + il> u = 4> u , Q+^o = 0, Q+^ = (6.27) 

Spectrum of the qubit-fermion supersymmetric system can be depicted in the following 
diagram 



Q 



E 



E = 



d+ ® /+ 



In the above example we use convention that (pi are even and ipi are odd, i = 0, u>. 

In terms of the supersymmetric quantum mechanics, we have here pair of nontrivial 
zero modes, and therefore the full even and odd spectra are identical. The analog of Witten 
index vanishes i.e. A = «q' -0 ^ — = 0. This might be surprising because, in 

conventional supersymmetric quantum mechanics even and odd spectra coincide except the 
zero energy ground state. However this effect is not new and is present in SUSY quantum 
mechanical models with a periodic potential [27, 28] or with local and nonlocal potentials 
[29] as well as in the model of spin | particle in a rotating magnetic field and constant scalar 
potential [30]. Here in the qubit-fermion system, the effect is of algebraical origin and 
comes from the structure of the model, not from the particular properties of the potential. 
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(c) qubit-qubit 

For completeness we rewrite in the same notational convention the Q- invariant two 
qubit system (considered already above). The Q-charges this time are even, but still nilpo- 
tent 

Qqq = iV^ d<g>d + , Q+ Q = -i^fuj d + (g>d (6.28) 
and graded commutation relations now read as 

[Qqq, Qqq}+ = " ( d+d ® 1 - 2N d) + (1 - 2N d ) ® /+/) = ffg> (6.29) 

h ( q F q = Qqq + QQQ' H QQ = H QQ + kH QQ (63 ^ 
[Hqq, Qqq}- = *H Q % [H QQ , Q+ Q ]_ = kH™ (6.31) 

(d) qubit-boson 

Qqb =iVub®d + , Q+ F = -iy/ujb + (g>d (6.32) 
and commutation relations now read as 

[Qqb, Qqb}- = w (-b+b® (1 - 2N d ) + 1 ® d+d) = #<J£ (6.33) 

H q F b=Qqb + Qqb, Hqp = Hq F + \Hq2 (6.34) 
For this system and A are even entities 

[H QB , Qqb}- = -XH§1, [H QB , Q+ B ]- = XH Q % (6.35) 
7. Separability and entanglement in Af-Hilbert space representation 

The commuting nilpotent variables were already used in the description of the entangle- 
ment in papers of Mandiliara et all. [10, 1 1]. Using advanced in the above sections notions 
of: ^-functions, ^-differential calculus and the ^-formalism of A^-Hilbert spaces we give 
a general characterization of pure state entanglement in terms of such a variables, and we 
find states considered in the Ref. [10] related to tanglemeters as a special case. 

As we have already seen, natural bases in tensor products of composite system are de- 
scribed by elementary functions, which encode automatically the "combinatorial" content 
of such non-simple (non-decomposable) tensors. In the approach presented in [10] there 
are discussed tanglemeters which are well defined for ry-functions with nonvanishing body. 
In addition we shall discuss the case of the //-function states with trivial body as well. In 
such a case one can apply the entanglement monotones introduced by Meyer and Wallach 
[23, 31] and their generalizations [32]. 

As it is known [33] one can also consider entanglement of fermions, bosons. In the 
case of fermions there is no Schmidt decomposition, but an analog of Schmidt rank (Slater 
rank) can be defined to classify entanglement in bipartite fermionic systems [34]. 
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(a) Factorization of rj-functions and pure state entanglement 

In this section we will show, how questions concerning the entanglement can be an- 
swered using the criteria of factorization of ^-functions. We will take inspiration from 
the factorization theory of functions of real variable, but of course proofs of the present 
theorems come from different reasoning then in conventional case. To begin consider the 
simple fact of linear independence of functions of one //-variable. This can be interpreted 
as special case of factorization aF{rj) + bG(rj) = => F(v) = — o^O. Functions 

F(rf) and G(rf) are linearly dependent iff the Wronskian of the following matrix 

W = ( FM ° M ) (7 1) 

1, dF( v ) dG{n) ) ' (/A> 

w = FqGi — F\Go vanishes. In particular functions with nonvanishing body, i.e. F, b ^ 

Ex 

can be written as F(r\) = F a e F o n . 

Let us define the entity H which will be important in expressing Wronskians for n-qubit 
systems. Namely, 

[§] 

U = EE(-!) fe (di k F(rJ)d In _ k F(rf)) | r7=0 (7.2) 

k=0 I k 

where di F(ff) — d$F{rf) — F(ff). In particular 

• for n = 2: 

H = F F 12 - F,F 2 , (7.3) 

• for n — 3: 

TL = F0F123 — F1F23 — F2F13 — F3F12 , (7.4) 

• for n = 4: we get the well known Cayley determinant 

H = F0F1234 — F1F234 — F2Fi 3 4 — F 3 Fi24 — F4F123+ Fi 2 F 34 + F13F24 + F14F23 (7.5) 

We do not distinguish in the notation the H for variuos values of the n, but it will always 
be clear from the context which one is under consideration. 



(b) Factorization of F(r]i,r] 2 ) 

As it is known [35, 36, 37, 16] a sufficiently smooth function of the real variables 
f(x, y) of the form f(x, y) — h(x)g(y) has to satisfy the d'Alembert condition i.e. 



d 2 lnf 
dxdy 



0. 



(7.6) 



This equation can be written in the form 



/ 

df 



df 



(7.7) 



But the set of solutions of (7.7) is larger then factorizable functions [36, 38] and there are 
various generalization of above equation for functions of several variables [37, 16]. 
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For the ^-functions analog of the relation (7.7) is also more general then that of the 
d'Alembert equation (7.6) moreover it is also necessary and sufficient condition for fac- 
torability of function. In case when F(r)i,r) 2 ) takes only invertible values, there exists 
/((^l) %)) = lnF(rn,r/ 2 ) and above conditions are equivalent. 

Let W12 denotes the following ry- Wronskian with respect to rji and r\ 2 variables 



Wi2{F{rn,^ 2 )) = detWu = 



F g£ 

dm 

OF d 2 F 



F diF 
d 2 F d 12 F 



— FqFi 2 — F\F 2 



dri2 9-qidri2 

(7.8) 

Let F(r)i,i] 2 ) be arbitrary function, then wi 2 (F) — iff F(i]i,r) 2 ) = G(f]i)G(f] 2 ), for 
some G and G. (Proof is given in Appendix C). Note that the Wronskian for the function 
of two rj variables has numerical values (i.e. its soul vanishes and there is no explicit rj 
dependence) and w\ 2 = H(cf.Eq.7 .2). 

As an example consider the Werner state which is represented by the function 

i>w{m,m) = + m), 

and we have that w\ 2 (ipw) = — \- For GHZ state 

wegetwi 2 (ipGHz) = |. but now there exists Iti{^/2i}jghz) — V1V2 and d'Alembert equa- 
tion also shows that ipGHZ is not factorable. Generalization of the d'Alembert equation to 
the system of n qubits was used by [10, 1 1] as a criterion of bipartite entanglement, and by 
means of the InF there was defined tanglemeter, as an additive measure of entanglement. 

Let us note, that from the condition of vanishing Wronskian we can distinguish several 
types of factorization: e Vl e V2 , 7]ie Vi , ■q i e nj , where i ^ j; i, j = 1, 2. On the other hand the 
non-vanishing Wronskian gives two types of ^-functions: GHZ-like ( e ±VlV2 , 1 ± r\ie ±Tli , 
e ±mm -|- e »j.) an( j W-like (1 ± rji ± rj 2 , ±771 ± rj 2 ). This can be seen from the easy to 
proof fact that for any function F(r/i,r] 2 ) there exists factor /(771) (or f(r] 2 )) and function 
G(r]i, r\ 2 ) such that 

F{m,m) = f(m)G{m,m), (7-9) 

where f(rji) = foe 7 " 1 ' 1 ■ Taking above ^-function phases with respect to both variables 
and assuming that / = 1 and g = 1 we get that 

F(m,v2) = f(vi)g(v2)G( m ,V2) = e^+^ 2 G( Vl , m ) (7.10) 

and the wi 2 (F) = wi 2 (G). Hence, both functions have the same factorization properties. 
For example functions related in this way to iPghz have the following form F(r/i,ri 2 ) — 
1 + hm + 92V2 + (fi92 + l)?7i?72 and these related to ip w are: F(r?i, r/ 2 ) = m + m + 
{.fi+ 92)mm- 

(c) Factorization of F {rji ,r] 2 ,ri 3 ) 

In the case of three variables we have to consider a set of the Wronski matrices, for 
all distinct pairs of variables. From such a bipartite information one can determine the 
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level of non/factorability of F(r]i,r]2, r] 3 ). This time the Wronskians depend explicitly on 
^-variables 

w 12 {F)(r} 3 ) = w 12 (F\ r , 3 =o) + {H + 2F3F 1 2)m = w 1 2(F\ r , 3 = )+H 3 m(7-ll) 
wi 3 (F){ V2 ) = w 13 {F\ V2=0 ) + (H + 2F 2 F 13 ) m = w 13 {F\ m=0 ) + H 2 r)2 (7.12) 
w 23 {F)(r h ) = w 23 (F\ m=0 ) + (H + 2F 1 F 23 ) m = w 23 (F\ m=0 ) + HivW-13) 

where H is given by (7.4) according to Eq.(7.2). What is interesting, one can express all 
terms in the above expansions of 77-function determinants in terms of the invariants of 
Wronski matrices. Namely, we find that the H k (F) can be written in terms of traces 

H k (F) = trW i3 {F\ m = ) ■ trW i3 (d k F) - tr {W l3 {F\ nk=0 )W l3 (d k F)) (7.14) 



Before we state criterion allowing to detect the existence of factorization (separation) 
of variables let us observe the following 

Weak factorization: Let ^(771,772,773) be such, that there exist functions G(r]i,r] k ) and 
G(rjj,rik) and ^(771,772,773) = G(r]i,ri k )G(rij,r) k ); i,j,k = 1,2,3, are all different and 
fixed, then Wij(F)(rj k ) = 0. 

However, inverse is not true, for example for the 77-function ^(771 , 77,773) = 7717/2 + 771773 + 
??2»73. w ij(F)(Vk) = 0, but there is no weak decomposition F = (r]i,r] 3 )G(r]2,r] 3 ). The 
non-vanishing Wronskians excludes possibility of factorization. 

One can introduce effective criterion selective enough to indicate full factorization of 
two subsystems (i — k)(j), with true separation of dependence on one variable. 

Strong factorization (separability): Let F = F(r]i , 772, 773), there exist functions G(i]i,r]k) 
and G(r]j) such that F (771, 772, 773) = G(rji,rjk)G{rjj) iff the following conditions are sat- 
isfied: Wij(F)(r)k) = 0, w k j(F)(r)i) = 0, Wij(d k F) = 0, and w k j(diF) = (i,j,k = 
1,2,3, are all different and fixed). 

The proofs of above theorem we present in the Appendix D. When there exists the loga- 
rithm of the 77-function then, condition for separability of variables takes very simple form 
and such criterion was already proposed in [10]. Namely, for n — 3, the dependence on 
rjj variable is separable iff didjln(F) = = d k djln(F), i,j,k — 1,2,3. Moreover, as 
it is easy to see, the vanishing of above derivatives of the InF is equivalent to the simple 
condition: w ik (F) = = w jk (F) e.g. did 3 lnF (r]i , rj 2 , rj 3 ) = (^13(^2) = and 
H 2 = 0). In this case, the rest of the conditions present in the strong factorization criterion 
is fulfilled automatically. Such simple characterization of separability in terms of logarithm 
generalizes to higher n [10]. However, the set of functions with vanishing body b(F) = 0, 
and hence not having logarithm, is very large and essential. It grows with the value of n. 
All Werner-like states follow into this type. 

From the point of view of the description of physical systems only the strong factoriza- 
tion gives information of the existing true bipartite separation of the systems, but even ex- 
istence of such a factorization still gives room for an another portion of non-separability in- 
side a two-qubit subsystem, except the totaly separable function i.e. F=G(rii)G(r] 2 ) . . . G{r} n ). 
The decomposition of F(ff) into factors is not unique in general. 
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Let us note that factorization properties of the 77-function F(t]i, 772, and its dual 
-kF(rii, 772, 773) are closely related, because 

w l0 {F\ m=o ) = Wij(dk(*F)), Hi(F) = H^F) (7.15) 

To illustrate how non/factorization properties of the F(rji , 772 , 773) are encoded in values 
of considered Wronskians let us consider some examples. For separability the normaliza- 
tion of the function is unimportant, but in view of further applications we shall consider 
normalized 77-functions. 
Examples: 

El: The Werner state is represented by the 77-function tpw = ^75 (^i + V2 + ^3), it is 
nonseparable, but it is non-factorable in the weak sense as well. We have 

Wii{^w\r, k =o) = ~, Wi = 0, (7.16) 
w t3 {d k ip w ) = (7.17) 

The cluster Werner state is represented by the 77-function dual to ij) W i.e. -kipw = 
(V1V2 + 77i77 3 + 772773), as before it is nonseparable and non-factorable in the weak 
sense. Namely, according to our observation concerning dual functions, 

Wij(*tp w \ Vk =o) = 0, 7Y 4 = 0, (7.18) 
Wij{d k ip w ) = (7.19) 

However, we should not be surprised to find that *ipw = ^ {4>w) 2 - 

E2: The GHZ state: tpcnz = ^75 (1 + Vi^rfo)- This function also is not factorable, even 
in a weak sense, because 

Wij(il>GHz\ Vk =o) = 0, ^k = \ (7.20) 
w^id^GHz) = 0, (7.21) 

Let us note that ipGHZ = *?Pghz, and nontrivial contribution to the Wronskian 
?7-function comes here only from the H i.e. Wij (tpGHz) = f\k- 

E3: In above examples functions cannot be expressed as a product of other functions 
except one factor is trivial, now let us illustrate case when decomposition into prod- 
uct of functions is not the factorization in a weak sense. Namely, let t/j = i(l + 

7?i7? 2 )(l + »?ir?3)(l + V2V3)) = \{ l + V1V2 + V1V3 + V2V3) 

WiM Vh=0 ) = \, H k = (7.22) 
Wij(d k il>) = ~ (7.23) 
and even no weak factorization exists. 
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E4: When we drop the first factor in above expression and normalize new functions, then 

-0' = ^5(1 + mm + mm) and 

Wl2(-0'U 3 =o) = = W 13 (ll>'\ m = ), W23(-0'| m =o) = ^, 

Hi = 

w 12 (d 3 ip') = -i w 13 {d 2 ip') = = ^(fliV'). ( 7 - 24 ) 

what here indicates weak factorization in 771 and 772 variables. One can express ip' 
as the product of two n = 2 GHZ states with common ^-variable . Namely, ip' — 

^(1 + + mm)- 

E5: The example of the state ip — (rj 2 + m + mmm) shows that assumption of the 
vanishing Wij (dkip) is important for the strong factorization. We have 

^12(^7,3=0) = = Wi 3 (V>|r, 3 =o), w 23 (ip\ m=0 ) = ~, Hi = 
wi 2 (d 3 ip) = ^ = wi 3 (d 2 tp), w 23 (d 1 ip) = 0, (7.25) 

Indeed, vanishing w\ 2 (ip) and wi 3 (ip) indicate presence of weak factorizations: ip = 
:^(l+f7i»72) (772+773) and ip = -^{l+mm)(m +m) respectively, but w 12 (d 3 ip) = 
w\ 3 (d 2 ip) + 1 show that the true separation of variables is not possible. Here we 
have product of rigged GHZ and Werner states. 

E6: Finally let ipbca factorable state (but not totaly factorable) e. g. ip = ^(1 + r] 3 + 

mm + mmm) 

w 12 (d 3 iP) = 1 U 3 = l - (7.26) 

wis^) = 0, 7i 2 = (7.27) 
w 23 (d 1 iP)=0, Wi = (7.28) 

In this case w i3 (ip) and u> i3 (djip) vanish, so dependence on 773 can be factorized and 
indeed ip <~ (1 + t]im ) (1 + 773 ) is a product of the GHZ-function and e m . 

Our classification of strong factorization of the functions for the n = 3 reveals the already 
known onion structure in the space of states discussed in [39, 40]. Firstly, we can distin- 
guish three sets Bi of mutually bipartite separable functions with nonempty common part 
of totaly separable functions. One can say that the Bi form a rosette which is surrounded 
by the set B\y of nonseparable functions of W-type nonseparability (cf. above examples) 
and of the set Bqhz of nonseparable functions of GHZ-type. It is interesting, that if we 
count dimensions of our 77-function spaces as follows: we take normalized functions; when 
function is separable we normalize each factor independently, then dimensions of above 
sets are the following 

dim(Bj) =3 + 1 = 4 (7.29) 
dim(n JJ B J ) = 1 + 1 + 1 = 3 (7.30) 
6xm{B w ) = 2 3 - 1 - 1 = 6 (7.31) 
dim(B GHZ ) = 2 3 - 1 = 7 (7.32) 



^12(^7,3=0) = 4, 
^13(^7,2=0) = 0, 
w 23 (ip\ Vl =o) = 0, 
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where for the Byv we have taken into account that there vanishes body F of the W-type 
function. These dimensions agree with the results obtained from the invariants theory [40]. 

(d) Factorization ofF{rji , 7/2, ??3, 

As before we shall consider Wronski 2x2 matrices to detect the possibility of weak 
factorization. The determinants Wij(F(rji, 772, rj 3 , 774)) are functions of two 77-variables 

Wij {F) = Wij(F\ Vk=ril=0 ) +H k {F\ m =o)Vk +ni{F\ m=0 ) m +Hki(F)vkm, (7-33) 

where e.g. H 3A = U + 2F 3 F 124 + 2F 4 F 123 - 2F 13 F 24 - 2F U F 23 and H is the Cayley 
determinant given by Eq. (7.5). We use notation 7Tj (F) — F\ Vi . Guided by the criterion for 
n = 3 we can state the following separability condition for n — 4: 

Strong factorization of 77-functions for n=4 (cases (i-j)(k-l) and (i-j-k)(l)): 



(i-j)(k-l) When the function F(rji, rj 2 ,rj 3 , 774) can be decomposed into the product of functions 
depending on separated pairs of variables F = G(i] il , r]i 2 )G(r]j 1 , r]j 2 ), I = i 2 } 
and J = {ji,j 2 } I U J = {1, 2, 3, 4}, where ii < i 2 , ji < j 2 then, 

Wi k 31 (F) = 0, w ik j, {d lk , F) = 0, w ikjl {d Jv F) = 0, w ikjl {d lk , d 3l , F) = 0, 

(7.34) 

where i k ^ i k , E I and jx ^ 3v € J- 

(i-j-k)(l) When a function F{r}\ , t]2, r/3, 774) has decomposition into the product F = G(rj il , rj i2 , t]i 3 )G{r}j ) 
then, 

Wi k j (F) = 0, w ikj (d lk , F) = 0, w ikj {d ik „ F) — 0, w ikj (d lk , lk „ F) = 0, 

(7.35) 

where J = {^1,^25*3}, J — {j} and as before IU J = {1,2,3,4}; i k 7^ i' k 7^ i k " ■ 
The proof is given in the Appendix E. For functions with nonvanishing body there exists 
logarithm ln(l+s(F)) and in such a case above conditions are equivalent to the statements 
that either 

d ik d jl ln{l + s{F)) = Q, k,l = 1,2 (7.36) 
for factorization (i — j)(k — I) or 

d lk djln(l + s(F)) = 0, k,= 1,2,3 (7.37) 

for factorizations (i—j—k) (I). We just get special cases of criterion considered in Ref. [10]. 
Let us stress that present criterions written directly in terms of Wronskians of ry-function 
F are more general and apply to all 77-functions. 

While for the n — 3 we do not have well defined 3x3 Wronski matrices, in the present 
case we can introduce nontrivial 4x4 ones in the following way 

_/ Wij (F) W l3 {d k F) \ 
^ \ WaidtF) Wij(did k F) J ' 
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where (ij) and (Ik) are ordered indices, hence i < j, I < k. Firstly, let us consider the body 
of these matrices, which can be written as conventional ones with the complex number's 
entries, in the following form 

It turns out that they are related to the known in the literature matrices obtained from the 
invariants theory, namely using results of Ref.([41]) we get 

delLyi = N, detL 14 , = M, detLf^ = L, (7.40) 

where Lf 3 T means partial transposition of the matrix L 13 

_ f W 13 (tt 2 tt 4 F) W 13 (ir 4 d 2 F) \ 

13 ^ w 13 (w 2 d 4 F) w 13 (d 2 d 4 F) J ■ 

The presence of partial transposition is important here, because invariants L, M, N satisfy 
the relation [41] 

L + N + M = (7.42) 

Such relation for Wronski matrices (7.38) would be troublesome. As it is easy to see, the 
77-function valued determinants of Cij give information about possibility of factorization 
and we are interested in all matrices dj, however for the complementary pairs of indices 
we get the same values of the determinant, so is enough to consider e.g.: L12, £13, L14. 
While above conditions of separability for n = 4 are invertible, we can get much weaker 
ones using determinants of 4 x 4 matrices: dj. It turns out that, a necessary condition for 
F(r}\ , 772 , 773 , 774) to be factorable into the product G(rj il , rj i2 , rji 3 )G(r]j ) is the following 

det£ ikj =0, k = 1,2,3 (7.43) 

When F is factorable into a product F = G(rj il , rn 2 )G(r]j 1 , ijj 2 ) then 

detCf? jt =0, M = l,2. (7.44) 

Let us note that for this type of separation of variables we take partial transposition of 
dj. To illustrate how above criterions differentiate various states let us consider some 
examples: 

El: For the n = 4 GHZ-state ipGHZ = ^(1 + J?i 772 773 774) we have 

w i: j(ipGHz) = ^ViVk, note that: b{w io {^ GHZ )) = (7.45) 
w ij {d^ G Hz)=0 (7.46) 
w iJ (d k di^GHz)=0 (7.47) 

and TL(tpGHz) = 5, but detCij = and also M = N = L = 0. The Wronskians 
Wij are sensitive enough to detect nonfactorability, and what is important, not only 
the body of Wij is important indeed. 
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E2: Similar result, we obtain for the n = 4 W-state \ (771 + 772+7/3+774). Here TL{ipw) = 
0, and 

Wij&w) = -l (7.48) 
«>y(fyM=0 (7.49) 
w iJ {d k d l il> w ) = 0. (7.50) 

This time, the nonzero content of the Wij is contained in the body of this Wronskian. 
As before: detdj = and M = N = L = 0. In conclusion we see, that here also 
the week criterion fails to detect essential nonseparability of the Werner state. 

E3: Present example shows that, the week criterion can be useful. Let us consider the 
following state: ip = ^(773 + 774 + 771772773 + 771772774). We have that detCf^ = = 
detC lk for i = 1, 2; k = 3, 4. But detC 12 = detC 3A = ^e^ 2 ^ 0. This signals 

actual factorability of this state, as ip = ^ghz^w^ m ^ remaining nonfactorability 
inside pairs of variables (1-2) and (3-4). 



8. Entanglement monotones and flavors of entanglement of 

77-functions 

The task of finding a good definition of proper, sensitive and operational entanglement 
monotone is not trivial. Despite that there are accepted monotones for 71 = 2,3,4 [42] 
there is still effort to modify them to suit special needs. Here we want to analyze known 
entanglement monotones in terms of 77-functions and relate them to the criteria of the 
non/factorization of these functions. What is remarkable, that the entanglement monotones 
for pure states can be expressed in terms of Wronskians of the 77-functions representing 
relevant state. Such representation of the entanglement monotones gives an insight onto 
the flavor of entanglement of particular family of states. 

In our discussion we shall follow the number of qubits and focus only on the pure state 
entanglement of bipartite systems. 

(a) n=2 

The well known entanglement monotone for this case, the concurrence can be ex- 
pressed using the defined above Wronskian 

^(771,772)) = 2^12(^(771, 772))!, <F,F>=1, (8.1) 

where the scalar product is defined in A/^-Hilbert space To{ff) and in components it takes 
the form < F, F >= |F | 2 + l^ol 2 + \Fi\ 2 + |F 2 | 2 + l^l 2 . Usingjhe notion of the 
comb [43] and antilinear mapping F F c , where F c = (a v ® a v )F one can express 
concurrence as 

C(F) = I < F C ,F > I = 2|F F 12 - FiF 2 |. (8.2) 
The operator a v ® a v is realized in the Fo(rj) in the following form 



-(did 2 + 771772 - mdi - r?i<9 2 ) 



(8.3) 
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It is interesting that by taking modulus of scalar product of the entries of defined above 
Wronski matrix i.e. 

V, = 2| < diF, F > | (8.4) 

and analogously defining 

V i = \<F,J i (F)>\, (8.5) 

where Ji(F (r/i, 772)) = F(-r]i,r]2) and J 2 {F(r] U r/ 2 )) = i^i, -772), we reproduce so 
called complementarity relations 

C 2 (F) + Vf(F) + Vf(F) =< F, F > 2 — 1, i = 1,2 (8.6) 

Using projections 7r fc | and 7Tfe|i we can express known in the literature parameters w, z, £ 
[44] by the ?7-scalar product and T>\ of a ^-function F 



w = V 1 (F(r ]1 ), F( m )) =J F( m ) A F( m )d m (8.7) 

F(m) A F( m )d m = V 1 {F{ m ) 1 F( m )) 
( = < F( m ), F(r h )>= J F^F^e^dfndm (8.8) 
z = < F( m ), F{ m )>= J F( V2 )F( m )e^dfj 2 d m (8.9) 
and we have thatC = 2|w|, Vi = 2\z\, V 2 = 2 j Ci - 

(b) n=3 

The 3-tangle [45, 46] is the principal entanglement monotone used to detect and mea- 
sure the "degree" of non-separability of 3-qubit systems. It makes use of the hyperdeter- 
minant known in invariants theory for a long time 

t 123 =4\Det(F)\ (8.10) 

It is sensitive enough to detect the GHZ state entanglement, which is maximal, but it ne- 
glects the entangled character of the Werner state. In the cited above paper, Coffman Kundu 
and Wooters derived relation for the 3-tange T123 and mutual concurrences of bipartite sys- 
tems of three qubits (we number qubits instead labeling them by the letter) in the following 
form 

^1(23) = ^12 + ^13 + T 123 (8-11) 

Using this equation, one can define entanglement monotone symmetric with respect qubit 
indices, averaged over possible configurations of qubits. Namely, 

I 2 2 2 2 

Q( F ) = 2^1(23) +^2(13) +^3(12)) = 3(^12 +^13 + ^23) + T 123 (8-12) 

Let us note that the Q is the n = 3 realization of global entanglement measure introduced 
by Meyer and Wallach [31] for arbitrary n (cf. also [47]for other form of this function). 
It can be seen that (cf. Appendix F for the proof) 

Det(F) = \ ]T(7^ _ 4^(^=0)^(9^), t 56 j + k (8.13) 
6 k 
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Using ^-realization of a matrices and the scalar product in M- Hilbert, one gets that 

< 4 j F,F>= -2{ Wij {F\ m=0 )+ Wij {d k F)) , (8.14) 

where bar denotes complex conjugation and ovf is a tensor product of I and a 2 matrices on 
i th and j th positions e.g. cr| 3 = I £g> 02 <E> cr 2 . Let us consider the following monotone 




(8.16) 



With the use of relation (7.14) the n(F) can be expressed solely in terms of determinants 
and traces of the Wronski matrices Wij (F), Wij (d k F) namely, 

= ^(2^|^(F| ))fc=0 )+ Wu (9 fc F)| 2 + |^((^(F| ))fe=0 ) (8.17) 

k k 

■trW l3 (d k F) - -tr(W l0 (F\ nk=o )W l3 (d k F))f - Aw l3 {F\ m=0 )w l3 (d k F)\). 

Above representation shows explicitly how entanglement and factorability properties are 
intertwined. 

From our previous discussion on Wronskians and Hi it follows that duality transforma- 
tion preserves the value of /j i.e. /J,(F) = /j,(*F). Let us test its behavior on some states. It 
will be instructive to come back to the functions considered in examples illustrating types 
of factorability of 77-functions in Sec c. Now the normalization of ^-function is important, 
to have fixed scale of values for fi. 
Examples: 

El: M^gfai +%+%)) = I 

The tpw state contributes to the value of fi only by the relative 2-qubit entanglement, 
the 3-tangle T123 vanishes. Analogously situation is for Mpw't ^(^(Vi^ + ViVs + 

mvs)) = § 

E2: The GHZ state cotributes to the value of fj, only through the 3-tangle ^{ipGHZ = 

^75(1 + mmm)) = 1 

E3: It is interesting to compare behavior of ji on states of the form ~ 1 + *ipw and 
— 1 + tpw For the first one we get £i(i(l + 771772 + 171773 + 172173)) = 1 with 
contribution solely from 3-tangle, but for the second state we get only contribution 
from the relative 2-qubit entanglement and /j(|(1 + 171 + 172 + 173)) = §■ 

E4: The function in this example and the one in the E5 give the same value of fi with 
almost the same mechanism, in both cases 3-tangle vanishes; //(-^=(1 + 771773 + 

%7? 3 )) = f • 

E5 = M-75(% + 773 + 771773773)) = § 
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E6: Here we have nonvanishing H3 and 3-tangle vanishes due to some cancelations of 
terms, again nonzero contribution to /j comes from the 2-qubit subsystem entangle- 
ment; = i(l + 773 + r/i?72 + V1V2V3)) = I- 



Note: Taking iffy = -j=(r)\+ 772) and ip G[ j z = ^75(1 + ^1^2) as the 3-qubit states with 
trivially factorized third qubit we get that n{ipyp) = § and ^(^q^z) = f ■ 



(c) «=4 



(i) Representatives of equivalence classes of entangled r]-wavefunctions 

Here we shall use results of the approach to entanglement based SLOCC classification. 
The problem of equivalence of the non-normalized pure states under permutations and the 
SLOCC transformations having determinant equal one was solved with use the invariants 
theory [42, 48, 41]. The representatives were found originally by [42] but we shall follow 
the modified form obtained in [48]. In terms of elementary 77-functions these nine invariants 
take the following form (for convenience of the reader natural form these functions is given 
in Appendix G) 



*i = 



a + d a—d. . . , b + c, . , 
— ^— e ' H — (cos(77i - 772) - cos(t7 3 + 774)) H — (cos(?7i - 773) (8.18) 

b - c 

- cos(t72 + 774)) H — (cos(?7i - 774) - cos(?72 + 773)) 

a + c-i - a-c + i b + c + i, 

— 2 — e h 2 — ( cos w i ~ w ~ C0S W3 + m)) h ^ — ( cos wi ~ 

b — c — 2 2 

- cos(?72 + 774)) H (cos(77i - 774) - cos(t7 2 + 773)) + -(sinfai + 7? 2 + 774) 

- sin(77 2 + 773 + 774) + + sin(77! + 773 + 774) - sin(7?i + 772 + 773)) (8.19) 
* 3 = ^e^ + w + ^-(cos( m - m) - cos( m + m )) + ^-(cos( m -774) 

- cos(?72 + 773)) + +i(sin(77i + m + %) - sin(77i + 772 + 774)) (8.20) 
*4 = ~Y- en + H COS (m - m) - C0S(772 + 77 4 )) + 7(cos(t72 - 773) - cos(t7i + 774) 

-2 sin 772 sin 773) + +^-^(cos(t7i - 772) - 003(773 + 774)) + ^(sin(?7i - 772 - 773 + 774) 
+ sin(?7i + 774) - sin(?72 + 773)) (8.21) 
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* 5 = - e " 1 " 2 +" 3r '*-27(sm77 2 cos(r /1 +r;3)+cos(r /1 +r;4)-l) (8.22) 
* 6 = i^V^* 1 ™ + a + l 2 +1 (cos( m - 773) - cos(t7 2 + m ) + a ~ l 2 X (cosfo - 773) - 

- cos(r7i + 774)) + -(sin(?7i + 773 + 774) - sin(??i + 772 + 774) + sin(77 2 + 773 + 774) - 

- sin(77i + 772 + 773)) + i(sin(77i + 772 + 773) - sin(?yi + 772 + 774)) (8.23) 

* 7 = sin 771 sin(77i + 773 + 774) + sin t? 2 sin(77i - 773 + 774) + i(sin(77 2 + 773 + 774) 

+ sin(77i - 772 - 774) sin(77i + 772 + 773)) (8.24) 

*8 = i(cos(«(77i + 772 + 773 + 774)) - sin(77i + 772 - 773 - 774) + sin(jyi + 772) - sin(77 3 + 774) + 

+ sin(77i - 773) + sin(?72 - 774)) + - (cos(t7i + 772 + 773 + 774) + cos(t7i + 773 + 774) - 

- cos(t7 2 - 773) - 2 sin(?7i + 7/3 - 774) + + sin^ - rj 2 + 773) - sin^ - 772 + 774)) (8.25) 
*9 = ie" 1+I ' 2 cos(7,3 - 77 4 )) + ^ e " 1 -'' 2 sin(7,3 + 774)) (8.26) 

where 77 = 774 = 771772773774. It turns out that <J>fc for k — 1, 2, 3, 6 are selfdual i.e. *(\E r fc) = 
^k- Let us observe that representation by means of the elementary 77 function is fluent in 
the sense that we have various identities for them, which can be used to change actual 
form of As noted by by Verstraete et all [49] the G a bcd state (here it is V&i function) 
is of generic type (any pure state of 4-qubits can be transformed into it) and with maximal 
4-partite entanglement. It is remarkable that using 77-functions it is easy to note that 

G abcd = tti = |e I ' 1 " 2 e I ' 3 " 4 + ^ e - r ' 1 " 2 e -*^ + ^(7 / i+772)(773 + 774) + ^(77i-772)(773 -77 4 ) 

(8.27) 

what means a sum of "diagonal" products of the GHZ and Werner states respectively of 
two-qubit subsystems of our total 4-qubit system. Namely, 

G W - *, - V 12) + -7^ (12) 7/> (34) + V 12 V 34) + -7/, (12 V 34) 

<^abcd — *i — 2^ghz+Vghz+ + 2^ ghz - Vghz - 2 w+ w+ 2 w w ~ 

(8.28) 

Above factorized form shows the flavor of entanglement of the *i function. 

(ii) Entanglement monotones based on polynomial invariants 

As it is well studied in the literature one of the main invariants in this case is the Cayley 
determinant. Another set of invariants consists of the determinants D uv = det(B uv ) [41]. 
Luque and Thibon give explicit form of them. It appears that in terms of the 77-functions 
they are related to the following Wronski matrices 

/ HfcnjF) H 3 {^F) H{-KidjF) \ 
Bij=\ Hi{KjF) Hij(F) MdjF) , (8.29) 
V HidiTTjF) Hj(diF) H{d t d 3 F) ) 

where the H, Hi and Hij are taken in appropriate form for 77-functions with n = 2,3,4 
variables. Taking determinants for pairs (ij) = (12), (13), (14) we get invariants which 
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are equal to D xy = det(B xy ), D xz = det(B xz ) and D xt = det(B xt ) in the notation of 
Ref.[41]. From the identity (7.42) it follows that 

L? + M 2 + N 2 = -2(MN + NL + ML). (8.30) 

To get symmetric (permutation invariant) monotones let us use the so called Schlafii basis 
{H, W, S, IT}, where W = D xy + D xz + D xt , E = L 2 + M 2 +N 2 and n = (L-M)(M- 
N) (N — L) . Osterloh and Siewert introduced for n = 4 the entanglement monotones \J-"i\; 
i = 1, 2, 3, 4, 5, which can be expressed in the Schlafii basis as follows (we modify slightly 
the form of relations from [43, 50]) 





= S(4W - 


-H 3 ) 




(8.31) 


T 2 


= 16(H 4 


-4HW- 


A(HD xt + ALM)) 


(8.32) 




= 32(H 6 


- 2<m 2 Y> 


- 6411) 


(8.33) 


Fa 


= 16(H 4 


-4HW- 


4(HD XZ +4LN)) 


(8.34) 




= 16(H 4 


-4HW- 


4{HD xy + 4MN)) 


(8.35) 



Therefore the for i = 2, 4, 5 are not symmetric, but can be combined into the sym- 
metric entanglement monotone. As proposed in [50] one can take a combination of them 
to get symmetric entanglement monotone 

= \T 2 + T A + F 6 \ = 16|3W 4 - imW + 8S| (8.36) 

Let us check how such defined symmetric entanglement monotone behaves on the follow- 
ing nontrivially entangled states: 

• the GHZ-state: ip GH z = ^(1 + VimmV^), *iPghz = iPghz 

• the Werner state: ipw = ^ (rji+ 1] 2 + m + V4:) and dual to it state *ipw = \( r limm + 

mmm + mmm + mmm) 

• the cluster Werner state: ip C w = 75(^1% + ViV3 + mm + mm + mm + mm)- 
It is a selfdual state *ipcw — ipcw 

. 0«)-family: 0<g = 4&wl + mm^ z± \ i.e. <l>% 4>™, <t>™, and <j>% 

. x W)-family: M = %^ z± + r?^Vf i), i-e. X^l X™. X ( ?L and x % 

. V«)-family: V>g = f (V&± + VkVi^z^l ™- ^ 

. A^-family: A<g = %^ z± + VkVl^z^ ^ ^+-> ^-+- 

. ^-family: $^ = f ( V , GHZ± +r ?j ^ ) +%^ ) ) e.g. <&^' 4) = f ty G HZ±- 

m^w 4) + m^w 4) ) 

• l> ± -states: $± = \($ghz± + VHipcw), 

where we take strictly ordered multi-indices (ij), or (ij, k) . Above families of states are 
spanned on pairs of the 2-qubit GHZ-states and W-states of the respective subsystems, or 
as in the case of $a±, <&± on total ipGHZ state and relevant ip^j), or ipcw states. From the 
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Table 3. Values of the entanglement monotones \T*2 \ and \Tz \ on selected famUies of states. 
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ry-function point of view they are natural and interesting as vectors in the .AAHilbert space. 

Some examples of them already appeared in the literature in the context of entanglement, 

(a) 

but without reference to the decomposition we use e.g. three states belonging to the (j)\ f_ 

family <^++\ 0++ ) were discussed in Ref. [50] and the ^2±^ was usec ' m tne 

Ref. [51] It is fruitful to compare behavior of j^l an d 1-^3 1 on above sets of states. We 
collect their values in the Table 3, together with the values of invariants contributing to 
the 1 an d 1-^3 1- Taking into account the way the final value is achieved one finds, that 

(ik I) ~ 

the cluster Werner state, -family and $± -family of states are interesting, because 

they give contributions from different channels: H, W, X, IT. The family $^ fc ^ is excep- 
tional, in the sense that, despite nontrivial contributions from all invariants: H,W 7 T, IT 
both monotones give zero as the entanglement measure. Other states are detected, if at all, 
only by one of available invariants contributing to entanglement monotones. The entangle- 
ment of the Werner state (or dual Werner state *ipw) is not detected at all. As we already 
discussed for the ipw> m that case one can get contributions directly from the Wronskian 
Wij(ipw) = — 4- ^ is worth noting that \!F' 2 \ an d I -^3 1 gi ye similar information about en- 
tanglement of considered states. Simultaneous vanishing of their values for arbitrary state 
is not obvious from the definitions of these entanglement monotones, except a particular 
case when the Cayley determinant H for a state vanishes and one of the determinants M, 
N or L is zero. Let, for example M^O and L = 0, then N = —M, and we have that 

|^| = 2 8 |M| 2 , (8.37) 

\T 3 \ = 2 12 |M| 3 , (8.38) 

or equivalently 

\F 3 \=2 4 \M\\^\. (8.39) 

We meet such a case for the families of states: V±^' ^±T' ^±±' X±±> as we ^ as f° r tne 
states dM k ' l \ 



9. Conclusions 

Description of the entanglement is a complex task, therefore any formulation which allows 
to see the problem from a new perspective is welcomed and may help in getting desired 
answer. The analysis of the qubit systems presented in this paper makes use of the new tool 
suitable for description of two level boson-like object - the nilpotent commuting variables. 
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In spite of the nilpotency, the property which such variables share with the Grassmannian 
ones (used in supersymmetric theories), the nilpotent commuting variables yield theory 
which is different then one known from supersymmetry. There, fermionic systems are de- 
scribed by means of the anticommuting, hence nilpotent variables with the use of the well 
developed supermathematics. Here, the differential calculus is not the supercalculus, but 
the new one, with modified Leibniz rule. For qubits, we do not want to use the anticom- 
mutation relations suitable for fermions, but qubit commutation relations which are, on the 
other hand, not bosonic ones, but define parafermions. Qubit has two dimensional repre- 
sentations and commuting nilpotent variables automatically provide correct properties of 
the single qubit and multiqubit systems. In analogy to the (super)symmetric theories there 
exists the classical mechanics for systems described by commuting nilpotent variables [5]. 
It is worth noting that, nilpotent variables emerge also in quantum field theory of complex 
systems. Supporting the idea that the nilpotent commuting variables give natural language 
for description of two level systems, but composed one - not fermionic. 

In the present paper there is developed further the formalism which can be used for the 
description of qubits. This formalism seems to be effective, practical and easy to interpret. 
Firstly, in this context nilpotent commuting variables were used by Mandilara at al. [10, 
11]. On the other hand, the present author developed the formalism suitable to describe 
classical systems with nilpotent commuting coordinates, in analogy to psedomechanical 
systems described by anticommuting variables. The peculiarities of appropriate differential 
calculus and an analog of the variational calculus were studied [5]. Hence, by now, various 
building blocks of the theory involving nilpotent commuting variables are at hand. In this 
work we have introduced the A^-Hilbert space of 77-wave functions allowing to formulate 
the Schrodinger description of quantum system including the analog of the Schrodinger 
equation (cf. [10]). It is remarkable that simple two qubit system, when not degenerate has 
the eigenvectors system consisting of Bloch vectors, and hence entangled. We have also 
addressed the question of hybrid, composed systems in analogy to the supersymmetric 
boson-fermion systems. As specially interesting we find the qubit-fermion system. It is 
really supersymmetric, in its simplest form it is composed of two, two level systems, having 
superalgebra of charges and supersymmetry transformations. It exhibits the effect present 
also in the conventional quantum supersymmetric systems with a periodic potential, that 
the Witten index vanishes. 

Many characteristic state vectors considered in literature in the context of entanglement 
have very simple and natural representation as ^-functions. They are just exponents, or 
trigonometric functions. In the present paper there are also introduced and studied the sym- 
metric ^-polynomials, Hermite ^-polynomials etc. The strong indication that //-functions 
are well tailored for representing the qubit systems comes also from the description of 
entanglement. Invariants known in the literature for various n = 2,3,4 brought to the 
physics from the classical invariants theory, frequently with complicated origins, here ap- 
pear as simple expressions in Wronskians of ry-functions. The questions of flavors of pure 
state entanglement and factorization within this formalism have natural, intuitive setting. 

We hope that the results of the Mandilara et al. [10] and results obtained in the present 
work show that the fy-formalism is suitable and effective tool for the description of en- 
tanglement. We think that it is something more then an ad hoc tool. The nilpotent classi- 
cal mechanics [5] and nilpotent quantum mechanics as well as field theory with nilpotent 
commuting variables [12, 9] reveal reach structure and deserve further study, although they 
describe non-fundamental particles, but composite ones. 
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Appendix A. Permanent and Hafnian 

As it is known there are four important functions from the set of square matrices into the 
ring or field over which matrices are considered. Namely determinant and permanent, then 
Pfaffian and Hafnian [52]. The permanent, like determinant is defined for any square ma- 
trix, Pfaffian for antisymmetric matrices and the Hafnian for symmetric ones, but in both 
cases n — 2k. The determinant and Pfaffian are sensitive to the parity of permutation while 
permanent and Hafnian are not. This makes great difference in properties of above objects 
and the level of complication in computing them. The determinant/Pfaffian is very univer- 
sal object and more easy to handle then permanent/Hafnian, its computational complexity 
is of polynomial type while according to the Valiant's conjecture: there is no polynomial- 
sized formula for the permanent [53]. This is important obstruction in calculating plane 
partitions in any models using bipartite planar graphs. Such a questions are extremely im- 
portant in e.g. dimer problem in statistical mechanics. 

Let A be an n x n matrix, the permanent of A = (a^ ) is a sum over permutations S n 

per(A)= II^W' (A1) 

<t€5„ 

so det is just a signed sum of the same type. Analogous relation is for Pfaffian and Hafnian. 
Let A 2nX 2n beasymmetric matrix and abe apartition i.e. ji), i 2 , J2), ■ • ■ , {i n ,jn)} 

(where ik < jk, h < h < • • • < "in) and a be a corresponding permutation 

_ / 1 2 3 4 ... 2n - 1 2n\ 

\ h h *2 h ■■■ ^ jn ) 

a ► A a = cn 1 j 1 di 2 j 2 ■ ■ ■ a i n j n (A 2) 

By Pfaffian one understands the signed sum 

Pf(A)= (~l) lal A a , (A3) 

q£S 2 „ 

and by Hafnian, Hf(A), 

Hf(A) = A «, ( A 4 ) 

a£S 2n 

The Pfaffian satisfies, among others, two important identities [52] 

Pf{Af = det(A), Pf(B T AB) = det(B)Pf(A). (A5) 
Analogous relations, in general, are not valid for the permanent and Hafnian (cf. [52]) 

Appendix B. Index conversion table 

To aid readers coming from various notational conventions we collect conversions of in- 
dices in binary, decimal and 77 function expansion component notation. 
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n=4 

binary decimal //-function 



aoooo 


a 


F 


aoooi 


a\ 


Fi 


aooio 


a 2 


F 3 


aoon 


«3 


F 3i 


aoioo 


ft4 


F 2 


Ooioi 


a 5 


F24 


oono 


a 6 


F23 


a 0111 


a 7 


F234 


OIOOO 


«8 


Fx 


OlOOl 


a 9 


F 14 


aioio 


aio 


F13 


CllOll 


an 


F134 


CfllOO 


ai2 


F12 


anoi 


013 


F X 24 


aino 


ai4 


F\23 


dim 


ai 5 


-Fl234 



Appendix C. Proof of factorization for n=2 

When F(r]i, 772 ) factorizes, then obviously wi 2 (F) — 0. Now, let F (771, rj 2 ) = F +Firjx + 
F2m + F X 2T]\r]2 and W\ 2 {F) = 0. 

• Let F , Fi, F2F12 7^ 0, then F 12 = ^A. Hence 

F(m,V2) = F (l + ^r/i + || 772 + ^p-??i772^ = F cxp j^r/i + ^772 



and, e.g. G( m ) = F e^ m and G(r7 2 ) = eW" 2 



(CI) 



• Let F a ^ 0, F i2 = 0, then e.g. F 2 = and Fi 7^ and F(t7i, 772) = 1 • (F + Fi?7i) 

• Let F = 0, then Ft or F 2 vanishes. Let F x 7^ 

F(vi,V2) = Ftf! +F 12 ViV2 = Vi(Fi + F 12 7] 2 ) (C2) 
Analogously for F 2 7^ 0; F = ?7 2 (F 2 + Fi 2 ?7i). 

Appendix D. Proof of the strong factorization criterion for n=3 

Let us consider the (12)(3)-separation. We have to show that there exisits decomposition 

^(771,772,773) = G(77i, 772)6(773) iff 

w 13 (F) = 0, w 13 (d2F)=0, (Dl) 
W2 3 (F) = 0, ™ 23 (<9iF)=0 (D2) 

Let F = G( m , 772)6(773) = GG 

G • G <9iG • G 
G-<9 3 G <9iG-<9 3 G 



wi 3 (F) 



(D3) 
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and 

d 2 G ■ G <9i<9 2 G • G 



wi 3 (d 2 F) 



d 2 G-d 3 G 9 1 d 2 G-d 3 G 



(D4) 



Analogously W2s(F) = = W2z(d\F). Now, let us see how factorization follows from 
vanishing Wronskians. We have 

w 13 (F\ m=0 ) = 0, ,H 2 = 0, w 13 (d 2 F)=0, (D5) 
W23(F\ m=0 ) = 0, ,Wi=0, w 13 {d 1 F)=0, (D6) 

(D7) 

Taking \(H\ + H 2 ) = and \(ft\ — H2) = we obtain the following set of conditions 

F F 13 -F 1 F 3 = 0, (D8) 

FoF 23 ~F 2 F 3 = 0, (D9) 

F F 123 -F 3 F 12 = 0, (D10) 

F 1 F 23 ~F 2 F 13 = 0, (Dll) 

F 1 F 123 - F 12 F 13 = 0, (D12) 

F 2 F\2 3 — F\ 2 F 23 = (D 13) 

(2) (2J 

For any F (771, 772, 773) we can write F (771,772, 773) =F (771, 773) + 772 F (771, 773) with 

(2) 

F (»h,%) = Fa + F1V1 + F3773 + ^13771% = F + F 1 r h + m (F 3 + F 13 i ll ) (D 14) 
and 

(2 J 

F (771,773) = ^2 + ^12771+^23773 + ^123771 773 = F2 + F12771+ 773(^33 + ^123771) (D15) 

To find the factorization of the F from conditions (D 9-D 13) we have to consider various 
forms of the F. The main families of solutions are labeled by F = (A) and F ^ (B) 
and obtained as follows: 
A: F ± . 

1. Fi,Fij,Fi 23 ^ 0. Using (D9) and (D 13) above relations (D 14) and (D 15) can be 
written as 

(2) F-l F 3 p p 

F= F e^ Vl ■ e^ m , F= F 2 e^ Vl • e^T" 3 (D16) 

£3. -F23 
and from (D 10) we get that, e F a 7,3 = e 7,3 . 

2. F 3 = 0. Then F 13 = F 23 = F 123 = 0. Hence, 

(2) Fi (2 - } F12 

F= F e 1! ^ m , , F=F 2 e^ Vl , (D17) 

therefore 773 dependence factorizes trivially. Let us mention here that above reason- 
ing in independent of the logarithm criterion, but (as discussed in the main text) is 
equivallent to the it. 
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B: F a = 0. 



(2) 



(2) 



1. F 3 = 0, Fi = 0; and F 2 ^ 0, F 13 = 0. Here F= and F has full expansion 
(D 15), but 773 dependence factorizes because of the Eq.(D 13). 



2. F 3 = 0, F = 0; and F 2 = = F 12 . Hence, 



(2) 



(2) 



F= Fi 3 ?7i?73, F= (F 23 + Fi2 3 ?7i)?? 3 
3. F 3 = 0, Fi = 0; and F 2 = 0, F 23 = = F13. Here we have 



(2) 

F=0, 



F= (F12 + Fi 23 77 3 )r7i 



(2) 



-Pi 3 



(D18) 



(D19) 



(2) 



4. F 3 = 0, F l ^ 0; and F 2 = = F 23 . In this case F= F ir]l e~ m and F = 

^"l23 -^13 

Fi 2 ?7ie~ p P~ T ' 3 . Due to (D 13) both functions have common factor e~ p T 7?3 . 

(2) 

5. F 3 = 0, Fi ^ 0; and F 2 7^ 0, F 13 = F 23 = F i23 = 0. Hence F= Fuji and 



F= F 2 + Fi 2 r/i, thus factorization is trivial. 



(2) 



6. F 3 = 0, Fi ^ 0; and F 2 , F 12 , F 13 , F 23 , Fi 23 ^ 0. We have F= F^e ^ 

(2) 



7J3 



then 



2 3 3 

due to Eq. (D 13) F = F 2 e~ F ^~ ??1 e^~ m . Both functions have common factor e - *^ 3 , 



because of Eq. (D 12). 



(2) 



7. F 3 7^ 0; and Fi = F 2 = F12 = 0. Here factorization is obvious, F= (F 3 



(2) 



Fi 3 ?7i)??3 and F= (F 23 + F 123 r]i)m 



Appendix E. Proof of the factorization condition of F for n=4 

Case (i-j)(k-l): Let F(r?i, r? 2 , %, = Gfai, Vj))G(r]k,Vi)- We have 



w 4 fc(F) = 



and also 



w ik (djF) 
w ik (djdiF) = 



GG diG-G 
G-d k G dtG-duG 



djG ■ G didjG ■ G 
djG-dkG d^G -d k G 







= 



djG ■ diG didjG ■ d t G 
dfi ■ d k d t G didjG ■ d k d t G 

Analogously relevant determinants vanish for other configurations of indices. 
Case (i-j-k)(l): Let F(r]i,r]2, 773, m) = G(Vi, Vj,Vk)G(r]i) then, 



(El) 

(E2) 
(E3) 



wu(F) 



GG d t G ■ G 
G-diG d l G-diG 







(E4) 
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and also, as above, 



Wa{djF) 



djG ■ G 



didjG ■ G 



djG-diG did.G-diG 



and 



wuidjdkF) 



djdkG ■ G didjdkG ■ G 

djd k G ■ dtG didjdkG ■ diG 

Analogous relations take place for other choices of indices i, j, k. 



(E5) 



(E6) 



Appendix F. n=3 Hyperdeterminant in terms of Wronskians 

The hyperdeterminant is frequently written in the form 

Det(F) = di - 2d 2 + 4d 3 , 



where 



di 
d 2 



^0^123 



F0F3F12F123 + F0F2F13F123 + F0F1F23F123 + F2F3F13F12 

+F 1 F 3 F 23 F 12 + F 1 F 2 F 23 F 13 
d 3 = FqF 23 Fi 3 Fi2 + F\F2F 3 Fi2 3 . 



Let us recall that 



Hi — FqFi2 3 + F\F 23 — -F2F13 — F 3 F\2 
H2 = FqFi 23 — -F1-F23 + -F2-F13 — F 3 F\2 
H 3 — FqFi2 3 — F\F 23 — F2F13 + F 3 F\2 



(Fl) 



(F2) 

(F3) 
(F4) 



(F5) 
(F6) 
(F7) 



and 



U>12(F\ n3= o) = FoF\2 — F1F2, 

U!i 3 (F\ n2=0 ) = F Fi 3 — FiF 3 , 
W23(F\ m= o) = F F 23 — F 2 F 3 , 



u>i2(d 3 F) = F 3 Fi2 3 — Fi 3 F 23 (F8) 
w 13 (d 2 F) = F 2 F 123 - F 12 F 23 (F9) 
w 23 (diF) = F1F123 - F 12 F 13 (F 10) 



Direct inspection shows that 



Y^iU = 3di - 2d 2 



and 



^2 w ij( F \r}k=o)Wij{dkF) =d,2- 3d 3 



(F ID 



(F12) 



fc=i 



Together with Eq. (F 1) we get relation (8.13). It is easy to get also analogous relation in 
the non-symmetric form e.g. 



Det(F) = H\- 4 U , 23 (F| J)1=0 ) W , 23 (d 1 F). 



(F13) 
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Appendix G. Representatives of entanglement SLOCC classes in 

naural form 

We present here the literal translation of states given by Chterental and Dokovic in [48]. 
Such representatives were firstly obtained by Verstraete et al. in [49], but in a slightly 
different form. 

a + d T a - d b + c b — c 

*i = — ^— e h — —(vim + mvi)^ — — (mm + mm) h — ^-(.mm + mm) (Gi) 

a + c-i -r a — c + i. b + c + i. 

*2 = — 2 — e + — 2 — (vim + mvi) -\ 2 — (vm + mm) 

b — C — 1 % 

h 2 — ( r?l7?4 + ^ 2?73 ) + 2*- 771 + f]4 + vimm + mmm -m-m (G2) 

-mmm - mmm) 

*3 = |e^ 1,)2+ ^ + ^(r /l r /3 + r, 2 7 74 ) + ^(r ?1 r /4 + r /2 r 7 3) + i(r / 3+r 7l 7 72 r/4 (G3) 

*4 = + fe^^ + 773774) + i{-mV3 + V1V4) + ^-^-(mm + mm) (G4) 

+^(»?2 + m + mmm + ^1%% -m-m- mmm - mmm) 
* 8 = ^ew+i™ - 2i{ m + f]im -mmm) (G5) 

* 6 = ^±l e ^+^ + i±l±!(„ 1% + r /2??4 ) + ^plfe^ + ^ 4 ) 
2 H - 1 2 — 1 z 

h — — (v3 + mmm) h — — (»?4 + mmm) - -^{m + m + mmm + mmm) (G6) 

*7 = mm + mm + mm- mm + mm + mmmm + Kmm + mm- m- mmm 

-mmm + mmm) (G7) 

*s = -^-^ -m- mmm) + ~if-(m + mmm - mm - mm) + + ??i??4 + »?i??3 
+mmm + m+ mm + mm + mmm) - Km + mm + mm + mmm) (G8) 

* 9 = Je^e" 2 " 3 " 4 + »?2 + r?3% + i(r? 3 + ??4 - r? 2 ??4 - t?2%)) (G9) 
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